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Abstract. During the past 15 years a new field has emerged, which combines 
superconductivity and spintronics, with the goal to pave a way for new types of 
devices for applications combining the virtues of both by offering the possibility 
of long-range spin-polarized supercurrents. Such supercurrents constitute a fruitful 
basis for the study of fundamental physics as they combine macroscopic quantum 
coherence with microscopic exchange interactions, spin selectivity, and spin transport. 
This report follows recent developments in the controlled creation of long-range 
equal-spin triplet super currents in ferromagnets and its contribution to spintronics. 
The mutual proximity-induced modification of order in superconductor-ferromagnet 
hybrid structures introduces in a natural way such evasive phenomena as triplet 
superconductivity, odd-frequency pairing, Eulde-Eerrell-Larkin-Ovchinnikov pairing, 
long-range equal-spin supercurrents, 7r-Josephson junctions, as well as long-range 
magnetic proximity effects. All these effects were rather exotic before 2000, when 
improvements in nanofabrication and materials control allowed for a new quality of 
hybrid structures. Guided by pioneering theoretical studies, experimental progress 
evolved rapidly, and since 2010 triplet supercurrents are routinely produced and 
observed. We have entered a new stage of studying new phases of matter previously 
out of our reach, and of merging the hitherto disparate fields of superconductivity and 
spintronics to a new research direction: super-spintronics. 


^ This is an author-created, un-copyedited version of an article published in Reports on Progress in 
Physies 78, 104501 (50 pp) (2015). lOP Publishing Ltd is not responsible for any errors or omissions 
in this version of the manuscript or any version derived from it. The Version of Record is available 
online at doi:10.1088/0034-4885/78/10/104501. 
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1. Introduction: historical background 

1.1. Superconductivity interacting with magnetism 

Quantum phenomena fascinate us and challenge our imagination for over 100 years since 
the theoretical foundations of quantum physics were laid. The two prime examples for 
macroscopic quantum phenomena are magnetism and superconductivity Soon after the 
discovery of superconductivity in 1911 in Leiden by Kamerlingh-Onnes mm it became 
clear that a magnetic field acts in peculiar ways on superconductors. Silsbee unified in 
1927 the observations of Kamerlingh-Onnes of a critical transport current density and 
a critical magnetic field |3], and in 1933, following a suggestion by von Laue, Meifiner 
and Ochsenfeld performed experiments showing that the new state of matter is a true 
thermodynamic phase and that it expels magnetic fields from its interior |1]. Following 
this pivotal discovery, Shubnikov pioneered in the 1930’s the field of interplay between 
superconductivity and magnetic field, culminating in the discovery of what is now 
called the intermediate state, or the Shubnikov phase [^ij^ where magnetic flux partially 
penetrates the superconductor in a well defined region of the temperature-field phase 
diagram. After the notion of a macroscopic quantum wave function was introduced for 
the superconducting state by Ginzburg and Landau in 1950 [^, Abrikosov theoretically 
explained in the 1950’s type II superconductivity, where a magnetic field penetrates the 
superconductor in a regular array of fiux fines carrying quantized fiux [Z]. Abrikosov’s 
vortex phase exists between the upper critical field Hc 2 and the lower critical field Hd 
in type II superconductors. 

Soon the microscopic theories of superconductivity by Bardeen, Cooper, and 
Schrieffer (BCS) in 1957 |8], by Bogoliubov in 1958 [9], and by Gor’kov in 1958-59 
mu brought all previous studies on a firm ground. 

Ginzburg, in an attempt to formulate a theory for ferromagnetic superconductors 
(he considered the possibility of superconductivity in gadolinium), concentrated in 
1956 on the electromagnetic (or so-called orbital) mechanism, where a suppression of 
superconductivity occurs via the interaction of the Cooper pairs with the vector potential 
of the magnetic field due to their charge m- In the orbital mechanism the magnetic 
field leads to an additional kinetic energy of the condensate, and when this energy 
exceeds the condensation energy, the superconducting state is destroyed. As shown in 
1966 by Werthamer, Helfand, and Hohenberg, the orbital critical field for conventional 
isotropic, diffusive superconductors in the absence of spin-orbit and spin paramagnetic 
effects is iiforb = ^c 2 |t=o ~ 0.7Tc[—dHc2/dT]T=Tc [E], where Tc is the superconducting 
transition temperature. 

An alternative mechanism, the exchange mechanism, was suggested in 1958 by 
Matthias, Snhl, and Corenzwit [13] in order to explain the variation of in lanthanum 
with rare earth impurities and the correlation between the appearance of ferromagnetism 

^ Unfortunately, Shubnikov’s work ended abruptly when he was arrested on the 6th of August and 
executed on the 10th of November 1937 during the Stalin-Yezhov terror on the basis of falsified charges 
of antirevolutionary activities; he was officially rehabilitated in 1957. 
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and superconductivity in ruthenides. They observed that it is not the dipole field of the 
effective moments of the rare earth elements that causes decrease of superconductivity 
in these systems, but the spin of the solute atoms; the does not correlate with van 
Vleck’s famous curve of /ieff dl, but rather with spin. In these cases an exchange 
interaction mediated by conduction electrons, responsible for ferromagnetism, occurs in 
a material that by itself is superconducting. The exchange interaction via conduction 
electrons tries to align spins in a ferromagnet, whereas the spins in a Cooper pair are 
opposite for the usual case of singlet superconductors. These antagonistic tendencies 
led to the so-called paramagnetic effect of pair breaking. 

In 1959, Anderson and Suhl showed in a seminal paper na that in a system with 
coexisting superconductivity and ferromagnetism the magnetism is modified due to the 
suppression of the spin susceptibility at small wavevectors. As a consequence, magnetism 
becomes under certain conditions inhomogeneous with a new, finite, wavevector. This 
new state was termed the cryptoferromagnetic state. The wavelength of the magnetic 
inhomogeneity is much smaller than the superconducting coherence length (the typical 
extend of a Cooper pair) and much larger than the lattice spacing, thus keeping 
neighboring spins almost parallel, however reducing the exchange field averaged over a 
typical Cooper pair volume. In real (anisotropic) materials the corresponding magnetic 
structure is often a one-dimensional magnetic domain structure. 

1.2. Pair breaking by paramagnetic impurities and external fields 

If one introduces paramagnetic impurities into a singlet superconductor, there is the 
effect of pair breaking due to scattering from paramagnetic impurities, which reduces Tc 
with increasing impurity concentration. Following previous works by Herring ng and 
by Suhl and Matthias na for the limit of small concentrations, Abrikosov and Gor’kov 
developed in 1960 a theory covering the full range of concentrations up to the critical 
value, when superconductivity is destroyed; the dependence vs. concentration is 
described in terms of a single “pair breaking parameter” p by the Abrikosov-Gor’kov 
formula m- In particular, these authors showed that the possibility of gapless 
superconductivity exists in metals with paramagnetic impurities (this peculiar state 
was further studied in 1964 by Skalski, Betbeder-Matibet, and Weiss PI, by Maki 
[20], and by de Gennes [21] )■ The Abrikosov-Gor’kov model, which employs the Born 
approximation, works well, e.g. for rare earth (except cerium) impurities. For transition 
metal impurities in superconductors Yu [22], Shiba [2S], and Rusinov [21] discovered 
within the framework of a full t-matrix treatment of the problem that local states (now 
called the Yu-Shiba-Rusinov states) are present within the BCS energy gap due to 
multiple scattering between conduction electrons and paramagnetic impurities. 

De Gennes and Tinkham noted in 1964 [25] that for a time-reversal invariant 
superconducting order parameter, pair breaking can be related to the asymptotic long¬ 
time behavior of the time-reversal correlation function, r]{t) = limj_j.oo(0^(O)0(t)), 
with the time reversal operator 0. An exponential decay, rjlt) = exp“^*/'’'®, with 
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tq some correlation time, corresponds to pair breaking with pair breaking parameter 
p = {2ttTtq)~^. If, on the other hand, a nonzero limiting value 0 < 77 < 1 appears, then 
this leads to an effective weakening of the pairing interaction. 

When applying a magnetic field to a superconductor, apart from the orbital effect of 
the penetrating field, there is also an appreciable Zeeman coupling between the electronic 
spins and the magnetic field. It leads to a splitting between the electronic spin bands 
similarly as in a weak ferromagnet the exchange energy leads to a band splitting. This 
effect has been experimentally studied in the early 1970’s by Meservey and Tedrow 
in a series of classical papers (see [2S] for a review). If the magnetic field exceeds a 
certain value, superconductivity becomes energetically unfavorable due to this Zeeman 
coupling. This limiting field is the Pauli paramagnetic limiting field or the so-called 
Chandrasekhar-Clogston limiting field, and was predicted independently in 1962 by 
Chandrasekhar iza and by Clogston [2H]- It amounts to = \/ 2 A(T — {))/gp,-Q, where 
A is the excitation gap in the superconductor at zero magnetic field, g the electron 
(yf-factor, and /tb the Bohr magneton. The ratio a = \/2Hor:h/Hp, where Horh is the 
upper critical field in the absence of the Pauli term, was discussed in 1966 by Maki, and 
is known as the Maki parameter ra- 

1.3. Ferromagnetic superconductors 

Coexisting singlet superconductivity and ferromagnetism is rare, and can be achieved 
either by finding suitable crystalline materials (classical examples are the rare-earth 
ternary compounds ErRh 4 B 4 [HO] and HoMoeSs [HI| in a narrow temperature region 
below the Curie temperature), or by introducing magnetic ions in a superconducting 
material that order ferromagnetically. In the latter case, at high impurity 
concentrations, the possibility to obtain coexistence between ferromagnetism and 
superconductivity was theoretically predicted in 1964 by Gor’kov and Rusinov [32] . 
Furthermore, for the case that ferromagnetism and superconductivity coexist or an 
external magnetic field is applied, Fulde and Maki [33] showed that the Abrikosov- 
Gor’kov formula holds with a modified pair breaking parameter. For an early review on 
magnetic superconductors see [34] . 

In addition to the above cases of conventional singlet superconductors, there 
has been discovered a large number of unconventional superconductors, which are 
distinguished from the conventional ones by the fact that they break additional 
symmetries, e.g. the lattice symmetry of the normal state, or the spin rotational 
symmetry. The latter is the case for superconductors exhibiting spin-triplet pairing 
(see figure . In particular ferromagnetic superconductors are under suspicion of such 
a pairing state, as equal-spin triplet pairing is not sensitive to the exchange mechanism 
in the way singlet pairing is. This refers to some heavy fermion compounds, like UGe 2 
[HS], URhGe [HS], UCoGe [HZ], and Ulr [HB], in which superconductivity can coexist with 
ferromagnetism, and which have been studied in the past 15 years p|. 

The main problem for coexistence between superconductivity and ferromagnetism 
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Figure 1. Two electrons, each of which has spin s = |, can combine their spin 
angular momenta Si and S 2 to build a pair with total spin S' = 0 or a pair with total 
spin S = 1. In the first case, the pair is in a spin-singlet state (shown on the left). 
In the second case, the pair is in one of the three possible spin-triplet states: with 
spin projection Sz = 0,±1 (shown on the right). The pairs with Sz = ±1 are called 
“equal-spin” pairs with respect to the spin quantization axis (here z-axis). The spin 
state \S,Sz) is fully characterized by the two quantum numbers S and Sz- We use a 
spin-vector representation to visualize pair angular momenta, where the expectation 
value of the cosine of the relative angle between the two spin vectors within a pair, 
(5,5,|si ■S2|5,5 ,)/v/^, is —1 for singlet states, and | for triplet states. Thus, 
angular momenta can be visualized as antiparallel for singlet pairs, whereas for triplet 
pairs they enclose in average a relative angle of ~ 70.(even for “equal spin” pairs). 


is that a strong exchange effect destroys superconductivity unless the pairing is of the 
triplet kind. However, not many superconductors support triplet pairing. Furthermore, 
such systems are typically p-wave superconductors BDI, which are very sensitive to 
pair breaking by normal impurities (conventional superconductors are insensitive to 
scattering from normal impurities, which is the content of a theorem by Abrikosov, 
Gor’kov, and Anderson [HI 112]) ■ For this reason new avenues have been chosen to study 
such systems for s-wave and d-wave superconductors, which are much more common in 
nature. These avenues combine the phenomena of proximity induced superconductivity. 
Cooper pairs with finite center of mass momentum, and odd-frequency s-wave spin- 
triplet pairs. 

1.4- Cooper pairs with finite center of mass momentum 

When the paramagnetic limiting field is smaller than the orbital critical field, then the 
possibility of an inhomogeneous superconducting state with finite pair momentum arises 
near Hp. Such a state was theoretically predicted in 1964 independently by Fulde and 
Ferrell [1H| and by Larkin and Ovchinnikov |11|. It is known in the western literature 
as FFLO state and in the eastern as LOFF state. 
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Figure 2. The FFLO effect: due to the exchange splitting of the energy bands (filled 
states are shaded blue in the picture) by ± J the Fermi surfaces split as well. A Cooper 
pair with opposite spins can be accommodated at the Fermi surface only on the cost 
of a finite center of mass momentum. Note that spin is quantized here in direction of 
the exchange field J, entering the Hamiltonian as "Hexch = — J • cr. 

Consider superconductivity in a weakly spin-polarized ferromagnetic material (or in 
an external magnetic field). Electronic spin bands are shifted in energy with respect to 
each other by an amount 2 J (which in general can be anisotropic): e{p) —>■ s{p)—J{p)-a. 
The exchange field can be related to an effective magnetic field via = J (for free 

electrons the magnetic moment is /r = pe < 0). We use the convention that spin is 
quantized in direction of the exchange field. The energy shift translates into a splitting 
of the Fermi surface for the two spin species (see Figure [^, which results e.g. for small 
J from the relation 

hvfQ = Vf{pf^-Pf^) = 2J{pf) ( 1 ) 

obtained by linearizing the dispersion relation around the Fermi energy [here the Fermi 
velocity Vf is assumed to be approximately equal for the two spin bands due to the 
small splitting, and Pf = {pj^-\-pf^/2]. Thus, the system only allows for opposite-spin 
Cooper pairs and built from electrons at the Fermi energy 

if they carry a finite center of mass momentum ±.hQ/2 = ^iPff — P/4,)/2 {hQ is the 
momentum of the Cooper pair). This can lead to an order parameter characterized by 
a linear combination of terms with different Cooper pair momenta 

Afflo(-R) = (2) 

V 

The preferred directions are chosen by the system (spontaneously or due to crystal 
anisotropy and boundary conditions). Fulde and Ferrell suggested an order parameter 
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characterized by only one wavevector Q, leading to a spatially homogeneous modulus 
and a spatially varying phase. If the order parameter is characterized by more than 
one wavevector Q^, then the FFLO state exhibits an inhomogeneous order showing 
a periodic structure. Larkin and Ovchinnikov proposed that the ground state may 
show a one-dimensional modulation (neglecting orbital effects). However, more complex 
structures are possible in the general case of anisotropic metals. In general, the pair 
amplitudes develop an [S' = 1, = 0) spin-triplet component, and the inhomogeneous 

modulus of the order parameter is accompanied by an oscillation of the magnetization 
around its average, determined by the product of the singlet and triplet pair amplitudes. 

Originally the state had been studied for a superconductor with ferromagnetically 
aligned impurities. However, it is known today that the FFLO state is very sensitive 
to disorder ga ig, and thus most likely is found in clean materials in an external 
field. Gruenberg and Gunther studied in 1966 the infiuence of orbital effects in type H 
superconductors on the possibility to observe an FFLO state and found that although 
orbital effects are detrimental to the formation of the FFLO state, such a state can 
exist at finite temperatures in the mixed state of a pure superconductor provided the 
Maki parameter a exceeds 1.8 HU. Interesting candidates are quasi-two-dimensional 
superconductors in a magnetic field applied nearly parallel to the conducting layers. 
The theory for such systems was developed by Bulaevskii gS] in 1973, by Burkhardt 
and Rainer [49] and by Shimahara gO] in 1994, and by Buzdin and Brison gl] in 1996. 
For quasi-one-dimensional systems theories in terms of soliton lattice solutions were 
developed by Buzdin and Tugushev in 1983 g2], by Buzdin and Polonskii in 1987 [53] . 
and by Dupuis in 1995 glj. In compounds with large spin susceptibility and high Hc 2 , 
as in heavy fermion and intermediate-valence systems, a generalized FFLO state was 
discussed by Tachiki and co-workers, where the order parameter is spatially modulated, 
and planar nodes of the order parameter are periodically aligned perpendicular to the 
vortices [55] . 

The FFLO state only appears below a certain temperature T*, which determines 
a tri-critical point in the (T, i7)-diagram at which the normal, BGS, and the 

FFLO phases meet. A generalized Ginzburg-Landau theory for the vicinity of this tri- 
critical point was derived by Buzdin and Kachkachi gg. A self-consistent calculation 
of the field versus temperature phase diagram and order parameter structures for the 
FFLO states of quasi-two-dimensional d-wave superconductors is presented in Ref. 157]. 
Recently, a symmetry classification of the pairing amplitudes in FFLO states coexisting 
with vortices was suggested gg. For reviews on this topic see Refs, gg EOl EB [S2|- 

The FFLO state has not been unambiguously found in bulk materials to 
date, although there are currently candidates under debate, as the heavy fermion 
material GeGoIns gSl EB ESI Eg, or some organic superconductors such as k- 
(bedt-ttf)2Cu(ncs) 2 gu EH Eg Eg. However, a similar state can be realized in 
superconductor-ferromagnet hybrid structures, where superconducting pairs leak into a 
ferromagnet due to the superconducting proximity effect, as we discuss in detail further 
below (for earlier reviews see gB EB EB EB EH]). 
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1.5. Proximity effect 

The pivotal role for realizing the FFLO state in such hybrid structures is played by 
the proximity effect between a superconducting and a normal conducting material. The 
theory for the proximity effect was developed largely by de Gennes [TSICTIEHIES], and 
by Werthamer m l8T] in the 1960’s. It describes the effect of penetration of Cooper 
pairs in a normal metal where the two conduction electrons (or holes) of the pair stay 
correlated with each other over a distance that depends on the amount of disorder in 
the normal material and on temperature. In a clean material, the penetration depth is 
determined by the thermal coherence length c = hvf /2'Kkffr. It gives the distance, a 
quasiparticle with velocity Vf travels during time h/2'KkffT.i before thermal decoherence 
sets in. For ballistic motion and atomically smooth interfaces, the penetration depth 
depends on the angle of impact of the quasiparticles, i.e. the Fermi velocity should be 
replaced by the Fermi velocity component in the normal metal perpendicular to the 
interface, Vf^^- The notion of a coherence length determined by the Fermi velocity and 
a characteristic superconducting energy scale was introduced by Pippard in 1953 |H2]- 

In the presence of impurities, scattering decreases the penetration of Cooper pairs 
into the normal metal. If the mean free path is £, and the motion is diffusive, then the 
quasiparticle travels in average between two scattering events a time r = (-/vf. For a 
random walk, the variance for the x-component is given by cr^ = where t 

is the elapsed time. Coherence is lost when the total path traversed during the random 
walk is equal to he. vff = ^n,c- If we define the diffusive coherence length by 

^N,d — obtain ^ = ^n,c^/ 3. Hence, the diffusive coherence length is given by 

^N,d = \/^n,c-^/ 3, or in other words, ^N,d = \/hD/2TTk-QT with the diffusion constant 

D = v/3g 

The proximity amplitude is proportional to the pair potential in the superconductor, 
A, and the product of the transmission amplitudes through the interface of the two 
particles comprising a pair. Inversely, the loss of Cooper pairs in the superconductor 
leads to the so-called inverse proximity effect of a weakening of the superconducting pair 
potential near the interface in the superconductor. 

1.6. Andreev reflection 

A very closely related phenomenon is the so-called Andreev reflection, discovered by 
Andreev in 1964 |H2] and by de Gennes and Saint-Jaimes in 1963/64 IHIIHS] (see Ref. 
[HS] for a review). It describes the correlations between particles and holes at the normal 
side of the interface due to penetration of pairs. In fact, Andreev reflection and proximity 
effect are two sides of one and the same coin. 

Let us assume that a spin 'f conduction electron in the normal metal with 
energy fi -\- e near the chemical potential /i moves in positive x-direction towards a 
superconductor. The excitation energy e is assumed to be within the energy gap A 

^ In d dimensions the factor of 3 is replaced by a factor of d. 





12 



Figure 3. Andreev reflection: An electron at crystal momentum pairs with another 
one at P 2 under the conditions that their energies and p^-component of momentum are 
opposite to each other. This leads to a finite Px component and a resulting center of 
mass momentum hQ^f2 of the Cooper pair. A hole is reflected, which is drawn in this 
plot at the negative momentum and energy of the missing electron (in the excitation 
picture all excitation states inside the Fermi surface have reversed momentum and 
energy and are considered as hole excitation with positive excitation energy). The 
resulting Cooper pair transfers its momentum partially to the interface between the 
normal metal and the superconductor (for the case that there is a Fermi surface 
or Fermi velocity mismatch), and partially to the condensate in the superconductor 
when entering it. The amount transferred to the condensate results in a supercurrent 
consistent with current conservation in the scattering process. 


of the superconductor. The electron can be transmitted into the superconductor only 
together with another electron with spin energy — and momentum approximately 
opposite, to build a Cooper pair that enters the condensate at the pair chemical potential 
2/i (after some momentum relaxation has taken place). For A <C /i it is possible to 
approximate the electronic dispersion E{p) ^ pj), with the Fermi momentum 

Pf and Fermi velocity Vf = v{pj). Within this approximation, the group velocity of the 
incoming electron is equal to Vf^ and its excitation energy s = Vf • {p — Pf) . 

Let us consider a Landau quasiparticle with momentum p^ = Pfi + Spi^ spin and 
excitation energy Si = Vfi • Spi. We will consider pairing with another quasiparticle 
with momentum P 2 = —Pfi + spin and excitation energy S 2 = —Vfi • 5p2 (see 
Figure]^. We require that S 2 = —Si> This leads to the relation Vfi • {Sp 2 — Sp^) = 0, 
which allows us to choose pj^^ such that Sp^ = Sp 2 = Sp hold^ The missing electron in 

^ Choose Pfi = Pfi + {6pi — Sp2)l2, then p^ — Pfi = P 2 ~^Pfi = = {^Pi + ^P2)/2- This is again 

(in leading order) a Fermi surface point due to Vfi • {p^i — p^i) = 0, i.e. the shift is tangential to the 
Fermi surface. 
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the spin-4- band is equivalent to a hole excitation characterized by energy = +£i, 
momentum ~ Pfi group velocity v{p 2 ) ~ —'^/i, and spin t (a missing spin 4-)- 
The process can be considered as a scattering of a particle with charge e, momentum 
Pf + 5p, and velocity Vf into a hole with charge —e, momentum Pf — 5p and velocity 
—Vf. The electron-hole pair gives rise to a current density of eVf + (—e)(—U/) = 2eVf, 
twice the current density of the incoming particle. The incoming electron (velocity Vfi) 
is approximately “retro-reflected” (velocity ~ — U/i) as a hole |83]. This process is called 
Andreev reflection, contrasting the familiar specular reflection. 

A possible small angle between the reflected hole and the incoming particle is due 
to variation of the Fermi velocity between the momenta p^ and Pf^ Sp, and is of order 
oi e/Ef, with Fermi energy Ef |H7]. Retroreflection is perfect for perpendicular impact 
if no supercurrent flows along the surface. In the Andreev particle-hole conversion 
process, energy and spin are conserved, however charge is not, and momentum is only 
approximately conserved. The charge 2e and (partially) the momentum hQ = 25p 
are transferred to the superconducting condensate: two electrons with opposite spin 
enter the superconductor to create a Cooper pair with non-zero pair momentum, which 
joins the condensate, leading to a (small) supercurrent for finite excitation energies. 
The momentum HQ is thereby partially transferred to the interface (provided a Fermi 
surface or Fermi velocity mismatch exists), and partially to the condensate resulting 
into a supercurrent consistent with total charge conservation. In equilibrium, and in 
the absence of a macroscopic supercurrent, this current is exactly canceled by the inverse 
process when a Cooper pair enters the normal metal and converts a hole into an electron. 

For an atomically clean interface, if there is a macroscopic supercurrent with 
Cooper pair momentum pg parallel to the interface present in the superconductor, then 
momentum conservation requires that the created Cooper pair enters the condensate 
with the required Cooper pair momentum parallel to the interface, such that 25p^^ = Ps- 
One obtains (we assume Vf^ > 0) 25px = {2s — Vf ■ pg)/vfx- In particular, for pg — Q 
only Px can change in the Andreev scattering process. In this case, as 5p^^ must be zero, 
the Cooper pair momentum is given by Qg = 2ejvfx (see Figure]^. 

The Andreev reflection process is phase coherent, i.e. the phases of the incoming 
particle and the Andreev-reflected hole are correlated. For a superconducting order 
parameter A = | A|e*^, the phase difference is <p{e)—x, between an incoming electron and 
retroreflected hole, and is ip{e)-\-x between an incoming hole and retroreflected electron, 
where ip{e) = — arccos(£/|A|) describes the phase delay due to the fact that during 
Andreev reflection the particle and hole amplitudes penetrate into the superconductor 
over a length scale = hv^p /\/|Ap — 

During their motion in the normal metal, the phase coherence is lost due to the 
slight difference in momentum of the electron and the hole, given by 25px = 2elvfx- 


^ For a supercurrent along the interface, x(i?) = x(i2o) + ^Ps ' (■P'~P'o), with li the spatial coordinate 
on the interface. The time the Cooper pair needs to enter the condensate is ^ = 0 V'^/ ^ which 
also determines the time delay between the maximum of an incoming electron wave packet and the 
maximum of the retroreflected hole wave packet, = hds^{e) = h/ \/| Ap — |88) . 
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When averaged over all momentum directions, the pair correlation function decays away 
from the interface algebraically ~ function of x, where we define = hvf/2e. 

For observables in equilibrium, the decay length can be obtained by replacing e in the 
expression for the momentum 5px by the lowest Matsubara energy IHH] ii^k-QT. This leads 
to an exponential decay on the length scale hvf /2TTkBT along each ballistic trajectory, 
which is precisely the coherence length from the proximity effect. 

For diffusive motion, non-magnetic scattering events conserve the time reversal 
symmetry, and thus scatter time-reversed states in a similar way. In particular a particle- 
hole pair with momenta both close to pj: will scatter in a particle hole pair with momenta 
both close to py (as scattering amplitude for pj. —)■ is the same as for —p^ —)• —p^ due 
to time reversal symmetry). Thus, scattering events change the momentum of particle 
and hole by the same amount, however do not destroy their coherence. This means, the 
dephasing between the particle and hole takes place on a semiclassical trajectory as in 
the clean case, however the trajectory changes direction multiple times in a random way. 
The diffusive process away from the interface is characterized by a diffusion equation 
for the pair amplitude /, given by dtf = Dd^f. With / ~ Qi(kxx- 2 etih) implies 
—2ie — —hDkl, resulting in a complex wave vector given for decay in positive x direction 
by kx = {i i: l)\/\e\/hD. This gives rise to an exponential decay on the length scale 
^/hDJ\^, accompanied by an oscillation on the same length scale. For observables in 
equilibrium we can replace e ii^k^T, leading to = —2'nk'BT/hD and an exponential 
decay on the length scale = ^/hD/27rkBT, again matching the length scale for the 
proximity effect in diffusive metals. The two phenomena, proximity effect and Andreev 
reflection are intertwined and cannot be discussed separately from each other. The above 
discussion assumes that the phase coherence is not weakened by additional effects. If 
phase coherence is weakened by inelastic processes, as for example magnetic scattering, 
this puts limitations on the proximity effect. In this case, one has to replace e by 
e = e -\- ias with the spin-flip scattering rate cTs = h/xg, where Tg is the corresponding 
life time. This means that in the ballistic limit exponential decay of pair correlations 
sets in on the length scale hvf/2as = u/Ts/ 2 (the factor 2 takes into account that both 
particles and holes are affected). The corresponding coherence lengths are obtained by 
replacing nk^T by xk^T -|- Og. 

A similar mechanism works for a weakly ferromagnetic normal metal, as for example 
a ferromagnetic alloy (see Figure]^ [HO]- In this case, the electronic bands are spin-split 
due to the exchange interaction by an amount of ±J <C Ef. The Andreev reflection 
mechanism now requires that the x-component of the momentum of the incoming 
electron with spin is pi^ = Pfx + (s^ + J)/vfx, which pairs with an electron with 
momentum p 2 x = —Pfx + (^ + This leads to a Cooper pair momentum of 

2{e^ J)/vfx. On the contrary, if the incoming electron has spin j,, it pairs with a spin t 
electron, and the sign of J reverses in the above expressions, leading to a center of mass 
momentum of 2(e — J)/vfx- Denoting Qj = 2J/hvfx and = 2elhvfx-, we see that 
instead of a singlet Cooper pair ('fj- — in the case of a non-magnetic normal 
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charge 2e 
momentum 

Cooper pair 

Qr2J/hvj^ Q^^2s/hVf^ 


Figure 4. Andreev reflection in a ferromagnet: an electron with energy c and spin t 
pairs with an electron with energy —c and spin resulting into a total momentum of 
h{Q^^Q j). For sufficiently small spin band splitting a Cooper pair of the form shown 
to the right results, which contains singlet and triplet components with amplitudes 
cos((5jx) and isin((5jx), respectively. The Cooper pair enters the condensate in the 
superconductor, transferring its momentum partially to the interface and partially to 
the singlet condensate, ensuring current conservation. The triplet component decays 
on the superconducting coherence length scale into the superconductor. 


metal, a Cooper pair of the form 


(t; 


o'i'Qjx 


-tte 


-iQjx^^^iQeX 


(3) 


is created. This is exactly the FFLO type of pair, as it was considered in the original 
work. Only it is now induced as a proximity amplitude instead of a bulk phase. It 
can be decomposed into a spin singlet (tJ. — 4,t) = |<S' = 0, = 0) and a spin triplet 

(ti + tt) = I'S' = 1, = 0) (where S is the total spin of the pair and Sz its projection 

on the . 2 -axis), leading to (see figure 

[cos((5ja:) |0,0) -|- i sm{Qjx) 11, 0)] (4) 


The proximity amplitudes in the ferromagnet are such mixtures. The triplet component 
also penetrates over a coherence length into the singlet superconductor, thus leading 
to singlet-triplet mixing in the supercondnctor close to the interface. The excess 
momentum during the Andreev reflection process is transferred partially to the interface 
and partially to the condensate (with the rate determined by the continnity equation 
for total charge conservation). 

In the diffusive limit the wave vector is given by Q± = (i + l)^y{e ± J)/HD (where a 
small positive imaginary part of £ determines the root), i.e. the pair amplitndes oscillate 
and decay on the same length scale. 
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singlet 


triplet 


singlet 


triplet 


singlet 




| 0 , 0 > 


/| 1 , 0 > 


i0,0> -/|1,0> 


| 0 , 0 > 


cos(^j^) |0,0> + i sm{Qjx) |1,0> 


Figure 5. Singlet-triplet mixtures between [S' = 0, 52 = 0) and \S = l,Sz =0) spin 
angular momentum states in an exchange field (or in an applied external magnetic 
field) in 2 ;-direction. Singlet and triplet spin states oscillate spatially out of phase 
relative to each other. Within the spin-vector representation the relative azimuthal 
angle between the two spin vectors within a pair varies as tt -|- 2Qjx, and its deviation 
from 0 (triplet) and tt (singlet) is a measure for the degree of singlet-triplet mixing. 
The top row shows a top-view of the bottom row. 


1.7. Josephson effect 


If two superconducting materials are brought in contact via a sufficiently narrow non¬ 
superconducting region (or via a weak fink), the macroscopic wave functions of the two 
superconductors may overlap. Under the condition that a phase difference exists between 
the superconductors, a supercurrent may flow between them directly through the non¬ 
superconducting region even under zero applied voltage. This effect was discovered by 
Josephson in 1962 IHH El and first observed by Anderson and Rowell in 1963 [93j . 
The Josephson effect is a hallmark of macroscopic quantum coherence, proving the 
macroscopic character of the pair wave function as originally predicted by Ginzburg 
and Landau |6]. Originally, the effect described the zero resistance tunneling of Cooper 
pairs between two superconductors through an insulating barrier, an SIS junction. The 
Josephson current from superconductor 2 to superconductor 1 (i.e. electrons are flowing 
from 1 to 2), I, is given by the set of Josephson relations 


I 

Ay 

dAx 

dt 


= F{Ax), F{Ax + 2 tt) = F(Ay), 
= X2-Xi - ^ J A-dr, 

= -^(^2- ^i) 

h 


(5 a) 
(56) 

(5c) 


where q* = 2e = —2|e|, xi c'-iid X 2 are the phases of the superconducting condensate 
wave function on either side of the contact, A is the electromagnetic vector potential, 
and q*$ is the pair electrochemical potential (<?2 ~ = U is the voltage across the 

junction). The function F describes the current-phase relation. The critical Josephson 
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currents in positive and negative flow directions are deflned via 
Ir 


ic+ = niaxF(Ax), 
Ax 


Ic- = — minF(Ax). 
Ax 


( 6 ) 


For non-magnetic barriers in the tunneling limit, Ic+ = Ic- = Ic and I = /cSin(Ax), 
where Ic is positive. The relations (5a)-(5c) are invariant against gauge transformations 
A' = A +VC, I'' = ^ — dtC, x' = X+ ^C- The function F{Ax) may not be single valued, 
for example for the case of a weak link multiple solutions can exist, and a hysteresis 
when sweeping Ay may appear. If a flux $ penetrates the junction, a characteristic 
Fraunhofer pattern appears when Ic is plotted vs. the flux. A multitude of interesting 
effects and devices based on the Josephson relations exist, and the reader is referred to 
review articles, e.g. [HiiinsiinniinziinB]. 

Here we are mainly concerned with the dc Josephson effect for the case when 
the non-superconducting material is an extended region of normal metal or metallic 
ferromagnet. With the help of the proximity effect via Andreev reflections, the 
Josephson effect can also occur in such devices. The current passing through the 
normal region is carried by the Andreev states, which build bound states below 
the superconducting gaps within the normal conducting region. The number and 
distribution of the bound states depend on details such as interface transmission, mean 
free path, and length of the normal metal. In general, there is a characteristic energy, 
the Thouless energy [21], given by hvf/L for the clean limit, and by HD/Ii^ for the 
diffusive limit. In normal metals coupled to conventional superconductors, there will 
be a low-energy gap in the spectrum of Andreev states, which for long junctions 
approximately scales with the Thouless energy and the transmission probabilities 
between the superconductors and the normal metal (possibly further reduced by inelastic 
scattering processes). This is the so-called minigap, found first by McMillan |l()()j . 
This minigap can be probed by scanning tunneling microscopy rnm \m\ . When a 
supercurrent flows across the junction, or when an external magnetic field is applied, 
the minigap is reduced and eventually closes. 

In the ballistic limit, each Andreev bound state with energy Fib.s. disperses as 
function of phase difference Ay, and the current is given by (g*/h)c?£'b.s./5^X- Apart 
from the current carried by the Andreev bound states, there is also a contribution from 
continuum states above the gap. 

The current-phase relation F(Ay) is in general strongly dependent on the details 
of the Josephson junction. Sinusoidal behavior is only observed under special 
circumstances, like for example in the tunneling limit. The critical Josephson current 
depends strongly on the length of the junction. For more details I refer the reader to 
the above mentioned review articles. 
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2. Symmetry classification of Cooper pairs in superconductors 


2.1. Consequences of Pauli principle and Fermi statistics 


Due to the Fermi statistics the pair correlation function, or the Gor’kov “anomalous 
Green function” ra. fulfills fundamental symmetries following from fermionic anti¬ 
commutation relations of the field operators. For fermions, the definition of the 
anomalous Green function |103i 1104] in Matsubara representation [89] is given by 

- F^p{ri,r2; r) = (T^ -iT)'l>^{r2, 0)}) 

= 6'(T)(^f„(ri,-ir)5';3(r2,0)} - 6'(-r)(^f^(r2, 0 )^f„(ri,-ir)} ( 7 ) 

where the function is defined as d(r) = 1 for r > 0 and 0 for r < 0, and where 
|r| < \h/kBT\. The function fulfills the fundamental identity 

FS^{r^,V2■,r)= -F^^{r2,v^■,-T). (8) 

which expresses the fermionic nature of the constituents of a pair. Going over to relative 
and center of mass coordinates, r = ri — r 2 , R = (ri + r 2 )/ 2 , the corresponding relation 
reads 


or after Fourier transformation in the relative coordinates r ^ p, r ^ z, 

F“j(p,z;R)= (10) 

Kubo-Martin-Schwinger (KMS) boundary conditions |l()5[n ()fi| allow for anti-periodicity 
in T with period h/k-QT'. T^(r, r — ih/k-BP] R) = —F^(r, r; R). This leads to discrete 
Matsubara energies z = Sn = (2n -|- l)7rA:BT, 

F“(p.£„;R)= -F«(-p.-£„;R). (11) 

We consider now spin-singlet (spin-antisymmetric) and spin-triplet (spin-symmetric) 
components: 

Ct = I =F F“) (12) 

(suppressing the spin indices hereafter), which fulfill 

F,^(p,£„;R)- F^{-p,-Sn;R) (13a) 

F,^(p,£„;R)= -F,’^(-p,-£„;R). (136) 


Each of these we can classify according to parity. We define even-parity and odd-parity 


functions 

F^{sn;R) = ^ (F^(p,en;R)±F^(-p,s^;R)) (14) 

(we suppress the momentum argument hereafter) and find 

Filt(£n;R)= K^(-Sn;R) (15a) 

Fil(£„;R)= -Fil(-£„;R) (156) 

F,^(£„;R) = -F,^(-£„;R) (15c) 

F,^(s^;R)= F,^(-s^;R). (15d) 
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Figure 6. Symmetry classification of Cooper pairs in inhomogeneous systems. The 
wavy lines symbolize odd-frequency symmetry of the pair amplitude. After [WlIT^ . 
With kind permission from Nature Publishing Group, and from Springer Science and 
Business Media. 

Thus, Fy(£„;R) and are even in Matsubara frequency, whereas Fy(£„;R) 

and F^{sn', R) are odd. 

After it was realized that in diffusive heterostructures with ferromagnets the type 
Fj^(£„;R) appears, which is according to the above a so-called odd-freqnency pair 
amplitude (for a review see [75]), a systematic classification according to symmetry 
of pairing correlations was undertaken along the lines explained above. Pair amplitudes 
are classified into four types according to their behavior with respect to frequency, 
momentnm (parity), and spin (see Fig. § 11071; 

Type A: spin singlet, even frequency, even parity 
Type B: spin singlet, odd frequency, odd parity 
Type C: spin triplet, even frequency, odd parity 
Type D: spin triplet, odd frequency, even parity. 

An identical classification scheme was independently proposed in Refs. [inHIIIDni- The 
fonr symmetry states above exhaust all possibilities compatible with Fermi statistics and 
the Pauli exclusion principle]^ The usual spin singlet, s-wave Bardeen-Cooper-Schrieffer 
(BCS) snpercondnctor [8] is of type A, while the spin triplet, p-wave superfiuid formed in 
^He [rnufTTn is of type C. Type D was first considered by Berezinskii ma in connection 
with early research on superfiuid ^He. Finally, type B was introduced in connection with 
nnconventional snperconductors by Balatsky, Abrahams and others |116[ I117[ I118| . 

Suggestions for realization of Type B superconductivity include generalized one¬ 
dimensional t — J and Hubbard models |119[ I12fij , two-channel Kondo models and Kondo 
lattices [Miinsinsiinii, and quantum critical points |I25|. Realizations suggested for 
the Type D state include disordered two-dimensional electron fluids in semiconductors 

® If more discrete degrees of freedom than the ones considered here are relevant (e.g. orbital) then the 
classification must be extended accordingly. 
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[miTTillTT^]. triangular antiferromagnets lEg, composite spin and orbital triplet 
superconductivity in two-channel Anderson lattices |I27|. Hund’s coupled pairing in 
double orbital Hubbard models imi, strong-coupling triplet superconductivity in 
Holstein Hubbard models ra, and one-dimensional systems with strong charge 
fluctuations m- 

2.2. Odd-frequency pairing amplitudes 

To date an odd-frequency superconductor has not been found in nature. However, 
in contrast to all cases discussed above, where the appearance of an order parameter 
due to spontaneous symmetry breaking (global phase symmetry in superconductors) 
was in the focus, in this review we are interested in the case of explicit symmetry 
breaking taking place locally in some spatial region, e.g. near interfaces, line defects, 
or inclusions. For this case an overwhelming number of experiments support the 
picture of the existence of proximity induced pairing states of type D. In its pure 
form (i.e. without additional components of different symmetry), it has been first 
predicted theoretically for a superconductor-ferromagnet proximity structure with a 
spiral inhomogeneous magnetization near the interface in Ref. |131) . Previous work 
had predicted the appearance of similar odd-frequency triplet anomalous functions in 
a model of coexistence of a superconducting phase and a helical ordering of localized 
spins, motivated by experiments on the superconductor ErRh 4 B 4 ra. 

In fact, the Pauli principle requires odd-frequency amplitudes to be present in any 
inhomogeneous superconducting state, not necessarily spin-polarized nm. For example, 
the case of a normal metal coupled to a superconductor involves pair amplitudes of 
Type B. Thus, odd-frequency amplitudes appear in heterostructures naturally due to 
breaking of translational symmetry at interfaces (similarly as Rashba spin-orbit coupling 
appears at interfaces in semiconductor heterostructures). They have been present in 
all treatments of inhomogeneous superconductivity for a long time, however were not 
explicitly named as such. Similarly, pair amplitudes of Type D appear as soon as 
spin rotational symmetry is broken, e.g. by an external magnetic field via the Zeeman 
coupling. If both spin rotational symmetry and parity are broken, e.g. at an interface 
between a superconductor and a ferromagnet, all four types of pair amplitudes are 
generated at the interface Similar conclusions have been reached in Refs. 

|108[ I109j . The important issue in discussing odd-frequency amplitudes in Ref. [131] is 
not the presence of Type D correlations per se, but their only presence. Thus, in diffusive 
structures those s-wave odd-frequency triplet pair amplitudes have a definite symmetry, 
and thus are of special interest for fundamental research. In regions where all amplitudes 
are mixed, simply all symmetries are absent, and a symmetry classification is not very 
useful. However, once such regions are coupled to reservoirs in which symmetries are 
asymptotically (far away from the spatial regions where the symmetry breaking takes 
place) re-established, then it is useful to discuss processes at interfaces in the light of 
those asymptotic symmetries. 
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To date it remains a great challenge to find direct experimental verification of 
the odd-frequency symmetry. There have been a number of proposals for indirect 
experimental verification via tunneling density of states studies, by inducing odd- 
frequency triplet superconducting correlations in a normal metal |134[ 113511136) . 

The odd-frequency pairing aspect in superconductor ferromagnet heterostructures 
was reviewed in Ref. |B|. A recent review dealing with odd-frequency pairing states 
and their relation to topological edge states can be found in Ref. |137j . Newer 
developments in finding odd-frequency order parameters include odd-frequency pairing 
states in multi-band superconductors |138l 1139) , at the surface of topological insulators 
|140) . in a controllable double-quantum dot system |141) . and in two-particle Bose- 
Einstein condensates [na, as well as odd-frequency density wave states [HHum]. 

Odd-frequency symmetry requires that the equal time pair correlator vanishes, i.e. 
electrons in a pair avoid each other in time. It has been argued that in such a case 
certain particle-hole symmetries based on a Green function approach are not appropriate 
and a Lagrangian formalism must be used |145[ 1146) . In particular, these arguments 
were brought forward for homogeneous odd-frequency states with an odd-frequency pair 
potential appearing due to spontaneous symmetry breaking. However, these arguments 
do not take into account that symmetries following from the full many body Hamiltonian 
of the system are fundamental. Should they be violated in a model Lagrangian approach, 
for which a corresponding model Hamiltonian cannot be found, then this simply means 
that the model Lagrangian is not appropriate for the problem in question m\- It 
has also been pointed out that spatially homogeneous odd-frequency states might be 
thermodynamically unstable in reality [Hi- 

3. Magnetically active interfaces 

3.1. Scattering phase delays 

Let us consider the basic quantum mechanical problem of a (quasi)-particle being 
reflected from an insulating region (which we assume at x > 0, potential V). Let us 
assume the particle has energy 0 < E < V with E = {k‘^ + k‘^^)/2m = V + + k'^^)/2m, 

and thus k{E) = [2m(y — E) + k(E) = [2mE — It is described by a wave 

function 4/(x,r||) = e**^i|i’ii ^-ikx^ at a: < 0 and 4/(x,r||) = at x > 0. 

The reflection amplitude for such a process is r = (fc — iK)/{k -\- in) = 
where k is the momentum component perpendicular to the interface and k determines 
the exponential decay of the wave function in the insulating region. In the limit for large 
n the phase of the reflected wave is flipped by tt with respect to the incoming wave. For 
finite K there is a phase delay with respect to this, given by <p{E) = tt — 2 arctan(fi:/A:). 
This phase delay results from the fact that the particle penetrates the insulator 
over a length scale h/n, and it corresponds to a time delay between the maximum 
of an incoming wave packet and the corresponding reflected wave packet |149) of 
Tci = hdEg^{E) = 2mh/(nk) (accompanied by a Goos-Hanchen shift along the interface 
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reflection amplitude 

{k-iK^I{k+iK^ = ^-2/arctan(«./^) 




reflection phase delay 
(p^ 71-2 arctan(A:yA:) 


spin rotation 


V 





Figure 7. When a Bloch wave is reflected from a ferromagnetic insulator, a phase 
delay appears between the reflected waves and the incoming waves, which differs for 
both spin directions. As a result, the spin-down reflected wave is delayed with respect 
to the spin-up reflected wave. If the incoming wave is polarized perpendicular to the 
ferromagnet’s magnetization, then the reflected wave has its spin rotated with respect 
to the spin of the incoming wave (see top part of the drawing). 


of s = 2hk\\/{Kk) |15()) ). 

3.2. Spin-mixing angle and spin-dependent scattering phase shifts 

We now assume that an exchange field J is present, leading to a contribution to the 
Hamiltonian given by 'Hexch = —Jo-. This can also be written as Hexch = —/t B^fi ■ o, 
with Seff = J/ 1^1 and /t is the (effective) magnetic moment of the charge carriers 
(negative for free electrons). Let us now consider a ferromagnetic insulating region, so 
that the two spin directions a have different reflection amplitudes — {k — iKa)l{k + 
iK(f) = e-2iarctan(Kcr//c)^ with ~ e in the ferromagnetic insulating region (see Figure 
[^. As throughout this review, we take the spin quantization axis along the exchange 
field, i.e. Kcr{E) = [2m{V — E — aJ) + such that the majority Fermi surface is 

assigned spin projection a = 10 The scattering phase delays now are spin dependent, 
= TT — 2 dxctan{K,^/k), and we can define a quantity •& = which is called 

spin-mixing angle imi or spin dependent interface scattering phase shift lEa in the 
literature (it is associated with a delay time = hdE'd{E) = 2mh/{n^k) — 2mh/{n^k) 
between the two spin components). For our example we have < Kp, and consequently 

® For free electrons (negative magnetic moment) this quantization axis would point opposite to the 
effective magnetic field. For a quantization axis in held direction, spin up and down would simply 
switch their roles. 
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Figure 8. Creation mechanisms for triplet pair correlations by direct and inverse 
proximity effect. See text for explanation. Reproduced and modified with permission 
from |153| . Copyright (2011) American Institute of Physics. 


> 0, > 0. Together with the spin-independent averaged phase ip = {p^ + ^^)/2, we 

write the reflection amplitude in a more general way as usa 

r{tp, t5, n) = -e**” ■ (16) 

with respect to a spin quantization axis h. The spin-mixing angle also describes the 
rotation angle of the spin components perpendicular to the axis h under reflection. This 
can be interpreted as a precession that the spins undergo around the axis h as a result 
of quantum mechanical penetration into the magnetic and insulating region. 

If one considers the Cooper instability in the presence of an interface with a 
ferromagnetic insulator, it becomes clear that pairing near the interface will be affected 
by these spin-mixing angles. In particular, a spin-up electron with momentum pointing 
towards the interface will have a relative phase shift with respect to a spin-down electron 
with momentum pointing away from the interface if the Cooper pair volume overlaps 
with the interface region. As the Cooper pair’s size is determined by the superconducting 
coherence length. Cooper pairs in a layer that extends a coherence length from the 
interface into the bulk will feel the spin-mixing phase shifts of the interface. This is 
shown in Figure]^ on the right. Near the interface, only a singlet-triplet mixed Cooper 
pair of the form 

(tT- T-t = cos(d)|0,0) -|- )isin(d)|l,0) 


(17) 
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can be present. According to this, we can consider the parameter d also as the parameter 
that governs the degree of singlet-triplet mixing at a spin-active interface. As seen in 
Figure there is a certain analogy between the creation of FFLO correlations in a 
ferromagnetic metal and the creation of singlet-triplet mixtures near a magnetically 
active interface in a superconductor. The role of the phase Qjx in the former is played 
by the spin-mixing angle d in the latter. As demonstrated in the middle column of Figure 
with increasing spin polarization of the ferromagnet in a superconductor-ferromagnet 
heterostructure, a shift from the creation of FFLO correlations in the ferromagnet to the 
creation of singlet-triplet mixtures in the superconductor takes place. This is because for 
very weak spin polarizations, the spin-mixing angle of the interface usually is of similarly 
small order and should be neglected in a consistent expansion in small parameters of 
the theory. For strong spin-polarization it is important and must be taken into account. 
The FFLO correlations in the ferromagnet show an inverse behavior: they are important 
for weak spin polarizations, however are restricted to atomically small distances from 
the interface for strongly spin polarized ferromagnets. 

The presence of spin-mixing angles has many consequences. For example, it leads to 
Andreev bound states at the interface, as predicted theoretically |154[ 11551 
m\, and verified experimentally |16n) . It also is responsible for giant thermoelectric 
effects in non-local setups HSU 1162] and is the main ingredient for creating triplet 
supercurrents in strongly spin-polarized ferromagnets |155L[T33] . It also crucially affects 
point contact spectra |lb3l 1164111581116511166111671116811169) . Finally, it is the main 
cause of the inverse proximity effect and the magnetic proximity effect in strongly spin 
polarized hybrid structures [EsiiiHsiinniiin]. 

3.3. Scattering matrix and induced triplet correlations 

The interface scattering matrix connects incoming Bloch waves with outgoing Bloch 
waves at an interface between two itinerant electron materials (metals, half metals, 
itinerant ferromagnets), or at a surface of such a material with an insulator. Relevant 
for transport are Bloch waves with energy close to the Fermi energy and momentum 
close to the Fermi momentum, in which case the Bloch waves are assumed to describe 
quasiparticle excitations. “Incoming” and “outgoing” refers to the projection of the 
quasiparticle’s Fermi velocity on the surface normal. We distinguish between electron- 
like and hole-like quasiparticles. For electron-like quasiparticles the projection of the 
group velocity on the Fermi momentum is positive, for hole-like quasiparticles it is 
negative. This means that electron-like and hole-like quasiparticles associated with the 
same Fermi momentum have opposite group velocity projections on the surface normal 
and consequently have different scattering matrices. For atomically ordered interfaces 
the crystal momentum component parallel to the interface, P||, is conserved, and then the 
relation between hole (h) and electron (e) like scattering matrices is <S'^(P|| ) = Se*(-P||). 

In relation to superconductivity, it is sufficient to consider the normal state 
scattering matrix for quasiparticles at Fermi momentum and Fermi energy. This results 
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from a systematic classification of all terms during a perturbation expansion in the small 
phase space volume associated with quasiparticles, which leads to quasiclassical theory 
of superconductivity |188[ I196j . In terms of the normal-state scattering matrix, 


S = 



(18) 


we can study the superconducting amplitudes induced by the presence of a singlet pair 
potential on either side of the interface. To simplify algebra and to gain some intuition 
we limit ourselves to the case of small pair amplitudes, in which case one can linearize the 
boundary conditions in the pair amplitudes. In our notation, Ri describes reflection into 
the superconductor, and R 2 describes reflection into the ferromagnet. All reflection and 


transmission parameters in Eq. (18) are 2x2 spin matrices. We consider the ballistic 
case, for which the reflected and transmitted pair amplitudes are given in 

terms of the incoming singlet pair amplitude = fo{i(Ty) by boundary conditions at 
the interface between the superconductor and the magnetic material, which in linear 
order read 


A. = A Ri {»,) ft; (19) 

/„“ = fc U (»„) h- (20) 


The simplest case is that of total reflection from an interface with a magnetic 
insulator. In this case we have (choosing the quantization axis appropriately, and with 
a scalar phase ip) 


Ri = 


0 


0 

e ^ 


The reflected amplitudes for this case follow from Eq. (19) and are given by 
foM = l/o cos(t?) + i/o sin(tf )<tJ ia„. 


( 21 ) 


( 22 ) 


This justifies the expression obtained by the simple arguments leading to Eq. 0- 
Next, we consider the case of arbitrarily large pair amplitudes in a ballistic 
superconductor with spatially constant pair potential. This is justified for sufficiently 
small spin-mixing angles t?. Incoming and reflected directions are parameterized by the 
polar angle 9p^, measured from the surface normal. The cosine of this angle is called /r. 
Defining even parity (symmetric) and odd parity (antisymmetric) functions with respect 
to the propagation direction, /“ = [/(/r) — /(—/i)]/2 and /* = [/(/i) -l- /(—/i)]/2, we 
obtain all four possible symmetry components as summarized in Table It can be seen 
that the corrections to the singlet amplitudes are ~ sin^('d/2). Thus, to linear order 
in 1 ? the pair potential stays unaffected. With increasing d the singlet pair potential is 
reduced, leading to corrections to the expressions in Table 

For the case of an interface with a strongly spin-polarized ferromagnet, one has 
to consider Fermi surface geometry, due to the presence of different Fermi surfaces for 
the two spin projections in the ferromagnet. Various cases can occur, depending on 
the scattering channel (parameterized in the ballistic case by P||, see figure 1^: for 
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symmetry 

even frequency 

odd frequency 

even parity 

Type A, singlet 0, 0) 

_ 7rOnA(l-sm2|) 

02—Apsin^^ 

Type D, triplet 1,0) 

Pg i7T£n^ ^ sin 7? 

■'to ~ a2|^|2sin2| 

odd parity 

Type C, triplet 1,0) 

— ZTrOnS^A 1 sint? 

Type B, singlet 0, 0) 

TrSnSA^Asin^l 

02—Apsin^^ 

02—Apsin^^ 


Table 1. Symmetry components of the interface amplitudes inside a ballistic 
superconductor in contact with a ferromagnetic insulator, assuming a constant singlet 
order parameter A. Here, = y^|A p + 5 ^, 'd = 'd(/i), and = sign(/i) {ji is the 
cosine of the impact angle). From Supplementary Material of Ref. |133| . 


example for reflection from the interface on the superconducting side, there can occur 
total reflection from the ferromagnet, total reflection of only one spin component from 
the ferromagnet (so-called half-metallic channels), or partial reflection for both spin 
components. Depending on the geometry of the Fermi surface in the superconductor, 
one, two, or all three cases may occur. 

For illustrative purposes we will consider a simple model of an interface potential 
of width d between a superconductor and a ferromagnet, with corresponding (free) 
electronic dispersions Vs + /2m in the superconductor, Vi — aJ + /2m in the 

interface region, and V^u — crJ + k‘^/2m in the ferromagnet. If the magnetization in the 
interface region is collinear with that in the ferromagnet, there will be no transmitted 
pair amplitudes, as for a strongly spin-polarized ferromagnet only pairs are 

supported, which both are not generated in such a system. There is, however, an 
induced triplet component on the superconducting side, with spin projection zero on 



Figure 9. Fermi surface geometry for a superconductor (left)-ferromagnet 

(right) interface. Various types of scattering events include total reflection (green), 
involvement of only one (blue) or of both spin bands (red). After |172| . Copyright 
(2009) by the American Physical Society. 
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Figure 10. Induced triplet components at a superconductor-ferromagnet interface 
with an interface barrier of width d. Model free-electron dispersions are: Vs + !2m 
(superconductor), Id — crJ + j2m (barrier), VpM — crJ + j2m (ferromagnet). The 
singlet amplitude in the superconductor is /o, the induced triplet component is ftz 
(and ftz/fo is purely imaginary). The Fermi wavevectors for the spin bands in the 
ferromagnet are fixed in (a) and (b). The various curves correspond to various Fermi 
wavevectors in the superconductor (as illustrated in the top left diagram in each panel). 
Wavevectors are normalized to ko = y/2mEF. For (a) Vi = 1.3, Vfm = 0-5, J = 0.2, 
Id = 0 (black dashed-dotted), 0.5 (green full line), 0.9 (magenta dashed); for (b) 
Vi = 1.4, VpM = 0.8, J = 0.3, Vs = 0 (black dashed-dotted), 0.6 (green full line), 0.9 
(magenta dashed); energies are in units of the Fermi energy Ep. (b) corresponds to a 
half-metallic ferromagnet, with only one spin projection itinerant (i.e. a Fermi surface 
exists only for a = 1). 

the spin quantization axis. It is shown in Figure 0 for the case of a strongly spin 
polarized ferromagnet in (a), and for the case of a half-metallic ferromagnet (with 
only spin projection cr = 1 itinerant) in (b). When the barrier is absent, there is 
no induced triplet component for the regions connecting wave vectors that support 
Bloch waves for both spin projections. Only wave vectors for which maximally one spin 
projection in the ferromagnet supports Bloch waves contribute to the induced triplet 
amplitudes in the superconductor. This shows that for high-transmission contacts 
half-metallic ferromagnets are most effective for including triplet correlations in the 
superconductor. If a spin-polarized barrier is present, triplet amplitudes are created also 
in the barrier region, and their sign depends on the relative size of the Fermi surfaces in 
the superconductor and in the ferromagnet. In the tunneling limit, the creation of triplet 
correlations happens mainly in the barrier region, such that they become insensitive to 
the Fermi surface geometry 

In order to allow for proximity induced triplet amplitudes in the ferromagnet, one 
needs an inhomogeneous non-collinear magnetization in the interface region. We study a 
simple model where the interface barrier is magnetized in direction perpendicular to the 
magnetization in the ferromagnet. We again assume a spin quantization axis in direction 
of the exchange held in the ferromagnet, which we take as the .s-axis, and a barrier 
exchange held pointing along the x-axis in spin space. In Figure pT] we show for two 
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Figure 11. Induced triplet components at a superconductor-ferromagnet interface 
with an interface barrier of width d, spin polarized in x-direction, perpendicular to 
the ferromagnet’s spin polarization in 2 ;-direction. Other parameters are as in Fig. 

(a), with Vs = 0, O.QF^i? for each thickness d. The components shown are: in 

the superconductor ff^/iifo) = (black); ff^/{ifo) = 

(orange); ffy /fo = (turquoise); in the ferromagnet /™/(i/o) (blue); 

(red). 


values of barrier thickness all triplet amplitudes generated in the superconductor and in 
the ferromagnet. For each value of barrier thickness we show three representative Fermi 
surface geometries. In the ferromagnet the spin bands support only triplet amplitudes 
of the form /™ and /™, whereas in the superconductor we show all three triplet 
amplitudes and 

It is interesting to note that for intermediate barrier thicknesses, also a component 
is generated, although the magnetization always lies in the x — z plane. This is due 
to the fact that spins polarized in z-direction, precessing around the ^-direction, develop 
a r/-component (and similarly for spins polarized in a:-direction and precessing around 
the z-direction). In a real gauge of the singlet pair potential, and are purely 
imaginary, and is purely real. As a consequence, only the former two contribute to 
the magnetization in the superconductor. 

Furthermore, depending on the barrier thickness, the vector of the triplet 
amplitudes in the superconductor changes direction rapidly as function of impact angle 
with respect to the surface normal. Only for tunneling barriers the direction of the 
triplet vector is dominated by the barrier exchange field. 

The equal-spin triplet amplitudes generated in the ferromagnet have opposite sign 
for the two spin projections. This reflects the fact that a triplet component with zero spin 
projection along the exchange held in the barrier region (a:-direction) decomposes into 
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equal-spin pairs with respect to the 2 ;-direction with equal magnitude and opposite sign: 
(t-i- + 4't)a; = (“ tt + Whereas in the tunneling limit the largest contributions 

to triplet amplitudes arise from near-normal impact, for thin barriers the dominating 
contributions arise from wavevectors that have for at least one spin projection evanescent 
solutions in the ferromagnet. 

For a general orientation of the interface barrier exchange field we parameterize its 
direction by polar and azimuthal angles, a and (/pj^ respectively, measured from the z- 
axis in spin space. The well-known transformation formulas for basis vectors quantized 
along the direction a, cp in terms of basis vectors quantized along the z-axis read 

ta,<^ = COS -e ~^2 + sin -e* 2 (23a) 

= - sin 2 + cos 2 4 ,^ (236) 

and can be used to find the transformation for pair amplitudes, for example 

in - n)a,v = in - u). ( 24 a) 

(n + n)a,^ = - sin a [e-*nn). - e*nu) J + COS a(n + ;t).. (246) 

It can be seen from this that if a triplet component of the form ( |246| ) is created in 
the barrier region, it gives rise to /™ and /™ correlations in the ferromagnet with 
relative phase of vr -|- 2ip. This leads to a non-trivial current-phase relation in Josephson 
devices, in which two singlet superconductors are connected by the two spin bands of the 
ferromagnet, each with their own separate current-phase relation m The azimuthal 
angle ip plays a crucial role in such devices [Ml [ISHl [US]. 

For the case of a long-wavelength magnetic inhomogeneity one has to deal with 
the spatial variation of the pair amplitudes, which can be done within the framework 
of quantum transport equations, which we discuss in the next section. In general, the 
scattering matrix should be calculated from first principles or from microscopic models, 
taking into account disorder, band anisotropy, and micromagnetics. 

4. Theoretical tools 

Theoretical treatments of superconductor-ferromagnet heterostructures can be 
separated into two types: microscopic treatments and quasiclassical treatments. There 
are two main methods, which both can be applied to the two types of theories: wave- 
function methods and Green function methods. Green function methods are based 
on an asymptotic expansion of the fundamental many-body Hamiltonian, whereas 
wave function methods work preferably with an effective mean-field Hamiltonian. 
Green function methods allow in an easier way for the treatment of disorder, and 
for generalization to strong coupling with bosonic excitations and to non-equilibrium 
situations. The resulting four theoretical frameworks with their corresponding equations 
of motion are shown in table The Bogoliubov-de Gennes equations [THl El lead 

^ In the following the scalar scattering phase ip will not appear anymore in this review, such that no 
confusion with the azimuthal angle ip should arise. 
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microscopic model 

quasiclassical 

wave function 

methods: 

Bogoliubov-de Gennes 
equations 

Andreev equations 

Green function 

methods: 

Gor’kov equations 

Eilenberger-Larkin- 
Ovchinnikov equations 


Table 2. Theories for inhomogeneous superconductivity can be classified into four 
types, as shown in this table. 


in quasiclassical approximation to the Andreev equations for the envelopes of the 
waves [HS]- Alternatively, one can start from the microscopic Nambu-Gor’kov matrix 
Green functions obeying the Gor’kov equations uniinii. These lead in quasiclassical 
approximation to the Eilenberger-Larkin-Ovchinnikov equations [ESIIITE], where the 
concepts of BCS pairing theory of superconductors [8] were merged with the concepts of 
Boltzmann transport equations within Landau’s Fermi-liquid theory cni. An extension 
to non-equilibrium was developed by Eliashberg |178j and by Larkin and Ovchinnikov 
jl79j . In quasiclassical approximation the Gor’kov equations result into envelope Green 
functions that vary on the superconducting coherence length scale and are free of 
irrelevant fine-scale structures on the Fermi wavelength scale. Dynamical phenomena 
are described in Green function methods within the Keldysh technique HHoi. 

In this review we concentrate on two, which have been mostly used so far in the 
literature: microscopic wave function methods based on the Bogoliubov-de Gennes 
equations, and quasiclassical Green function methods, based on Eilenberger-Larkin- 
Ovchinnikov equations and their counterpart for diffusive systems, the Usadel equations. 
The two other frameworks, not covered in this review, are formulated in terms of 
Andreev equations and in terms of Gor’kov equations. 


Jj..!. Bogoliubov-de Gennes equations 

The Hartree-Fock-Bogoliubov (HFB) mean-field Hamiltonian 

^HFB ~ d^rd^r' (Hra,r'a' V’LV’r'cr' 

<T<T'e{th} 




(25) 


is the basis for the Bogoliubov-de Gennes theory. Hermiticity requires iLro-.r'o-' = 

and Fermi statistics leads to Ara,r'a' = ~^r'(T',rcr- The Hamiltonian can be rewritten in 

a compact form using the 4x4 matrix 


H(r,r') = 


[Hrcr,r'a']2x2 [^ra,r'a']2x2 

[^r<T,r^(T^]2x2 [■^r(T,r^(T^]2x2 


(26) 
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where, for example, 

\_Hra,r'cr''\2x2 


Hriy'f Hriyi 


(27) 


With the dehnitions d''f(r) = ^r^), ^(r) = i'iprt,Ah'>Py,'ipliV (here T 

denotes a transpose) the Hartree-Fock-Bogoliubov Hamiltonian reads 


'Hhfb J J ^d'^(r) i7(r, r') ^(r')^ + const. 


(28) 


This Hamiltonian is diagonalized by a Bogoliubov-Valatin transformation |181l 1182) . To 
achieve this, we introduce the notation for the Bogoliubov amplitudes {a e {+,—}) 


Ur,, 


r = 


^ ^rt,nQ; ^ 


\ 


5 ^nQ:(^) — 

?;* 

Y\.^na 



«rt,na 

y ^r4,,nQ; j 


\ Kl,na / 


(29) 


Here, the index a may, e.g., refer to another spin basis, or to a helicity basis if strong 
spin-orbit interactions are present. The quantum numbers n, a fully characterize the 
eigenstates of the set of Bogoliubov-de Gennes equations given by 


dVi7(r,r')t/„«(r')= ^naG„„(r) 

dh'H{ry)Unair') = -F;„„G„„(r). 


(30 o) 
(306) 


For each eigenvalue Ena > 0 there exists an eigenvalue —Ena, with corresponding 
eigenvectors obtained by the substitution Uj,a,na ^ Ka,na^ Vj,a,na ^ '^ra,na- One can 
divide the quantum numbers n in disjunct sets Z, V and Ef (and relabel them if 
necessary) such that ioi n & Z we have Ena = 0, for n G P we have Ena > 0 and 
for n G W we have Ena < 0. The eigenvectors can be chosen to build an orthonormal 
set, e.g. for n,n' within V 


I ^ ^ ^ ^ (^r(T,naW(T,n'a' ^^cr,na'^ra,n'a'^ ^nn'^aa', 


(31o) 

(316) 


We build for n G T* a 4x4 matrix t/„(r) from the four orthogonal eigenvectors 
t/,+ (r),G„_(r),t7„+(r),t7,_(r) (ora 4x2 matrix if only one value for a belongs to n). 
Similarly, for n & Z we build a 4 x 2 matrix Un+ (r), Un- (r) (or a 4 x 1 matrix if only one 
value for a belongs to n). Then the completeness relation 

5^f/,.(r)f/t(r')=i(r-r')i (32) 

nev 

(with the 4x4 unit matrix i, and with P = Z UV) holds, which reads explicitly 

^ ^ {^ra^na'^r'a',na cr',na) ^ ) 

n^V^OL 


(33a) 









32 


^ ^ ('^rcr,na'^r'cr',na cr',ricy) — 0 (336) 

ne.V,a 

{sn = 0 for n ^ Sn = 1 else). With this, the Hartree-Fock-Bogoliubov Hamiltonian 
is diagonalized by 

= '^Un{r)%, %= f dhul{r)i>{r), (34) 

nev ^ 

with, e.g., 7 ^ = ( 7 ),+ , 7 y, 7 „+, 7 „_), 7 ^ = denoting Bogoliubov 

quasiparticle operators and En = dia,g{En+, E^-, —-E'n+, —En-) corresponding energies, 
leading to 


ff(r,r')= 5;7(r)47(r') 

neV 

(35) 

T-LuFB = 9 7n-^n7n- 

n&r 

(36) 

Explicitly, Eq. (|m 1) reads 


t^rcr = ^ ^ (Ura,na'lfna T ^n’^rajUa'^na) 

neV,a 

(37o) 

-fr,a = dhY^ iKa,na^ra + C,naV'L) • 

a 

-0 

r- 

CO 

The Majorana property 7 ^^ = 7na holds if Vj.a,na = 

'^rancd ioV 71 E Z 8i SolutioU of 


the Bogoliubov-de Gennes equations can always be chosen of this form provided this 
is compatible with the boundary conditions. The single-particle part of i^ro-.r'o-' 

has typically the form 


Hra,r'a' = {^(P) “ + g(P) ' “ h(r) • <t| 

J aa' 


(38) 


with P = —iES/r — eA(r). Here, £(P) is even in P and g(P) is odd in P. The last 
two terms in Eq. ( |^ describe spin-orbit coupling and Zeeman coupling, and is the 
electrochemical potential. The order parameter is given self-consistently by 


A 


r(7,rV' 


(fvi(fr\ V^cr'aia{ (r, r', ri, 


rifJiy 


(71 cr'^ 


(39) 


with 


(^r 


''^ra) — 2 


^ra,na ^r'a'^na 


yUr'a' 


n^V,a 


,na J" 


tanh 


En 


2>k^T 


(40) 


The spin-dependence of these equations simplifies when the pairing interaction can be 
split into singlet and triplet parts, 

and similarly for the order parameter, 


^ra^r'a' (^^2)(7(7' H“ ^r,r' ’ ctct' ' 


(42) 
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The gap equation for spin singlet and triplet ferromagnetic superconductors has been 
scrutinized by Powell, Annett, and Gyorffy 

In order to study odd-frequency amplitudes one uses a Heisenberg representation 

Aa{t) = («r<T,naeS^'*“*7na + SnC,nae“ (43) 

ne.V,a 

and evaluates the pair correlation function |184j 

= ii^raWr'a'iO)) (44) 

using (qLTna) = f{Ena), (bnaTna) = l-/(^na), with/(£■„„) = [l-tanh(£’„„/2A:BT)]/2 
the fermionic distribution function. Concentrating on local correlation functions, one 
obtains |184[ I185| 


(«rt,na'yr^,na TM4,na<t,na)Cna(^>^) (45o) 

n^V^OL 

Faair,t) = Fra,ra{t) = Y “r<T,na<a,naCa r) (45&) 

ne.V,a. 


for the singlet (F_) and the three triplet {F^, Fn) components, where 

CJi. T) = cos (^) tank (^) - i sin (^) . (46«) 

Ct.(i,T) ^ cos (^) - 1 - *sm (^) tank (^) . (466) 


The triplet correlations show a time dependence according to and consequently 
vanish identically for t = 0, whereas the singlet correlations, governed by 
survive. This expresses the odd-frequency nature of the local triplet pair correlations 
in contrast to the even-frequency nature of the local singlet pair correlations. In 
particular, for t —>■ —ir with h/k^T > r > 0, and applying KMS boundary conditions 
Cnai-^^ + i^/kBT,T) = the relations Q^{iT,T) = ±Q^{-iT,T) follow. 

Similarly, measurable quantities like magnetization, current density, and charge 
density can be expressed in terms of Bogoliubov quasiparticle operators using Eq. (37a) 
and evaluated using the solutions of the Bogoliubov-de Gennes equations. 


4-2. Quasiclassical theory of superconductivity 

The treatment of a superconductor in the presence of a Zeeman spin-splitting 
induced by an external magnetic held has been theoretically investigated in detail by 
Alexander, Orlando, Rainer, and Tedrow in 1985 M within quasiclassical theory 
of superconductivity. This theory, a generalization of previous studies for superhuid 
^He [wi rrmi \m\ . includes Fermi liquid effects self-consistently, as well as impurity 
scattering (including ballistic and diffusive limits), internal exchange helds, and spin- 
orbit effect^ It is formulated in terms of generalized Landau parameters 

^ The derivation in this work is more rigorous than in the later work by Dernier et al. in 1997 m- 
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as well as singlet and triplet pairing interactions based on Leggett’s 

treatment for clean systems [wmiD], and combines previous theories of Fulde dHHl 
and Buchholtz and Zwicknagl dso]. First experimental tests of these theories were 
performed by Tedrow and Meservey HU. Generalized Landau parameters lead for 
example to anisotropic renormalizations of the quasiparticle magnetic moment by 
coupling to electrons outside the phase space regions where quasiparticles live (“high- 
energy electrons”). 

For reviews of quasiclassical theory of superconductivity see Refs. [Ml IMl 
11881 11941 11951 11961 11971 11981 11991 12001 1201] . Quasiclassical theory is an expansion 
in quantities like temperature or gap divided by Fermi energy, or Fermi wavelength 
divided by superconducting coherence length [18811202112031120411205] . It predicts its 
own breakdown in low dimensions [ 206 ]. 

The central quantity in quasiclassical theory of superconductivity is the 
quasiclassical 4x4 matrix propagator or Green function g(pj:,R;e), which depends on 
the spatial coordinate R, on energy s, and on the momentum directions pj on the Fermi 
surface (out of equilibrium there is in addition a time dependence and the formalism 
can be extended to 8 x 8 Keldysh matrices m)- The 2 X 2 Nambu-Gor’kov matrix 
structure of g reflects the particle-hole degree of freedom mm- Its matrix elements 
are 2x2 spin matrices (if the spin-splitting of the energy bands is not comparable to the 
energy band width, otherwise they are scalar). We will use a notation of unit matrices 
and Pauli matrices in spin space and particle hole space, where (ctq, cri,< 72 ,(J 3 ) refers to 
spin and (tq, ti, T 2 , 73 ) to particle-hole degrees of freedom. We expand the elements in 
Nambu-Gor’kov space into spin-scalars and spin-vectors, 

gao + g-a ifcro + f-a-)ia 2 

i<y2{f<yo - /-cr) (^2{g(^Q -g-<T)a2 

with e.g. g = {gx,gy,gz) the vector part of ^n, and g its scalar part. Of particular 
interest for this review are /, /, describing spin-triplet correlations, and g, g describing 
the spin magnetization and the spin current that develop as a result of the generation 
of triplet correlations. 




4-3. Eilenberger equations 


We discuss first the case when the spin splitting of the energy bands is small (comparable 
to the superconducting energy scales), so that it can be treated as a perturbation around 
the un-split Fermi surface. In this case one can integrate out the energy dependence of 
the Green’s functions identically for the spins t and I m and assume spin-independent 
Fermi velocities Vf{pf). The propagator g obeys the Eilenberger transport equation 

[rTHlITTHI 


en -h, g 


ihvfVng = 0 


with the normalization condition nza 


(48) 


-2 2 - 
g = -TT To, 


(49) 
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where Vf{pf) is the Fermi velocity for the direction at the Fermi surface, and with 
h hgxt 3“ dmf -|- djmp (50) 


where hext are external potentials, are Fermi liquid mean helds (both diagonal and 
off-diagonal), and dimp is the impurity self energy. The matrix h(pf,R;e) entering the 


Eilenberger equation (48) has a similar structure as g, 


h 


( lyao + vf-cr (Actq-|-A- cr)i(T2 
zcr2(Acro - A-cr) a2{i>CFo - i>-a)a2 


(51) 


The Nambu-Gor’kov matrix structure contains some redundancy, which results into 
symmetries between the particle and hole elements |188) . They are expressed by the 
tilde operation q{Pf, R, s) = q{—Pf, R; —e*)* for q G {g, /, g, /, u, A, i/, A}. Matsubara 
propagators are obtained by (e ten) (where £„ = 'Kk-QT{2n + 1) is the Matsubara 
energy), retarded propagators by (e —)■ e -|- i5), and advanced propagators by (s —>■ 
e — iS). Further fundamental symmetry relations are: g(s*) = g(s)*, g(s*) = g(s)*, 
f(e-) = f(er, f(e‘) = -/(e) * (we omitted here for brevity the other arguments). 

The self-consistency equations for the impurity self energy dimp for impurity types 
i with impurity potential hj and impurity concentration n* is 

^imp(Pf, R;e) = '^ niii(pf,pf, R; e) (52) 


with the quasiclassical T-matrix equation 

= HPf,Pf) + 

+ {NfiPf)MPf,Pf)9{Pf^ R, £)ii{Pf,Pf, R, £))#.; (53) 

(.. .)p^ = J (fpj:/(2TThY ... denotes a Fermi-surface integral, and Nf{pf) = (27rd|'y/|)“^ 
is for a fixed Fermi surface point the (one-dimensional) density of states in perpendicular 
direction to the Fermi surface, per spin projection, in the normal state. 

The Fermi-liquid mean-field self energies d^f are self-consistently determined from 


i^mi{Pf,R) = k^TY, {Nf{Pf)W^{Pf,Pf)g{Pf,R-,e^))p> (54a) 

Sn 

UmiiPf^R) = k^TYiNfiPf) {pf,p'f)g{p},R-,e^))p, (546) 

A„,(pp R) = keT Y, (A'/(p;) V(P/, p',)f(p;, R, £„))„; (54c) 

£n 

A^t{Pf,R) = kBTYiNfip'f) {Pf,p'f)fm,R,e^))py (54d) 

£n 


Fermi-liquid interactions are parameterized by dimensionless Fermi-liquid parameters 

A"{Pf,Pf) = NfW^{Pf,Pf) and {Pf,Pf) = Nf 1T“ {Pf,Pf), with Nf = {Nf(pf))p^ the 
density of states per spin projection at the Fermi level, and the superconducting pair 

potentials by singlet and triplet pairing interactions and (where kF“, and V* 
are in general tensors). 
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Finally, Dgxt contains the coupling to external fields, for example the electro¬ 
magnetic coupling to a magnetic vector potential A{R), 

^orbital {Pf,R) =-eVf{pf)-A{R)T3 (55) 

(e = — |e|), and the Pauli-coupling to the quasiparticle spin 

Vspm{Pf,R) =-B{R)-{Pf)-S (56) 

with the effective quasiparticle magnetic momenf|^ 

^eff (Pf) = [i- {Pf)]l^e (57) 

where A® (pj:) = Nj' {Nf{p'j^) ^4“ {Pf,Pf))p'^, and the spin matrix in particle-hole space 


s=(^ ° ]■ (58) 

y 0 -Cr2<T(T2 J 

Pe is the electron magnetic moment (pe = ~|he 1)0 The vector potential A and hence 
the magnetic field {B = V x A) are calculated from the current density as 


V X V X A{R) = /ioj(^) 


(59) 


with the permeability of free space fiQ. The scalar electrochemical potential follows from 
local charge quasi-neutrality, and is zero in equilibrium. The charge current density is 
obtained from g{pf,R;en) via 

j{R) = 2ekBT {Nf{Pf)vf{pf)g{pf, R; £„))-^, (60) 

the spin current density for spin projection along the axis €„ via 

J^{R) = 2ek^T {Nf{pf)vf{pf)[ec, ■ g{pf,R-en)\)p^ (61) 

£n 

and the spin magnetization via 

M{R) =Xn B{R) + 2kBT Y (NfiPf) ^efr(P/)&(P/,-R;£n)}p^, (62) 

Sn 

where Xn— 2A^/(1— Aq)ij,1 is the normal state spin susceptibility, defined in terms of 

the parameter Aq= {Nf(pf)A^ (P/))py- The local density of states for a given spin 
direction e is calculated at real energies as 

NeiR;e) = -^lm{Nf{pf)[g{pf,R;e+) + e-g{pf,R;e+)])p^ (63) 

with £_|_ = £-|-zfi (5 > 0 infinitesimal). The Eilenberger equation for g, its normalization 
condition, the equations for the self energies a, and Maxwell’s equation for A, must be 
solved self consistently by iteration together with the appropriate boundary conditions 
imposed on the propagator and the vector potential |2fi7[ I2fi8[ 1209112in[ I211j . Boundary 
conditions in quasiclassical theory are notoriously difficult to derive, as the quasiclassical 


® i denotes a 3x3 unit tensor. 

Pe = —ffMB/2, with the Bohr magneton gs = \e\h/2meC and g = 2.0 
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Figure 12. In quasiclassical approximation boundary conditions involve relations 
between the envelope functions of Bloch waves on the two sides of an interface. These 
envelope functions in general show a jump, even when the Bloch wave functions 
themselves are continuous. 


propagators show in general jumps at interfaces (see figure [^, in contrast to microscopic 
propagators which are continuous. 

A powerful way to implement boundary conditions and to solve Eilenberger 
equations is the Riccati method [na EH EH EH. Modern versions of boundary 
conditions for Eilenberger equations using these techniques are presented e.g. in Refs. 
piillTHilETHlETH] . 

Instead of applying an external field B{R), a superconductor can be spin-polarized 
via the inverse (or magnetic) proximity effect when in contact with a ferromagnetic 
material. This effect already appears when the ferromagnet is insulating. A theory for 
the spin polarization and the associated induced exchange field in the superconductor 
has been developed by Tokuyasu, Sauls, and Rainer in 1988 HSH. Although differing 
in details, the general mechanism presented there is essentially the same as for 
all subsequent studies of the inverse proximity effect in superconductor/ferromagnet 
heterostructures in the last decade: the appearance of a Cooper pair spin polarization, 
or of triplet pair correlations, in the superconductor creates a finite spin magnetization 
inside the superconducting region, extending roughly a coherence length away from the 
interface. This decay length is dictated by the decay of triplet pair correlations as the 
bulk of the singlet superconductor is approached. 


4.4- Normalization condition and transport equation in spin-space 

The interrelation between triplet amplitudes and induced magnetic moment can be 
understood within quasiclassical theory already by the normalization condition for the 
propagator. The normalization condition (49) encodes important information about 
the spin structure of the diagonal and off-diagonal propagators izrzi. In equilibrium, 
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physical solutions require the condition 

Tr (g) = 0, i.e. Qq = -go (64) 

which also ensure local charge neutrality to leading order in the Fermi liquid expansion 
parameters. It is useful to introduce the notation g± = {g &)/2, i'± = {u ± i^)/2, 
u± = {u ± u)/2. The external and internal fields are distributed according to: fiorbitai 
in Uspin in which also contains the internal exchange field J] spin-orbit band 
splitting in In terms of measurable quantities, g_^_ determines the magnetization, 


g_ the spin current density, and go the charge current density. One obtains from (49) 


9o + ^" = ff-f-f-9l-9- 
2gog- =if X f 

‘2go9+ = If - ff 


with g\ = g^^ 


9y,± 


9z,± = 9± 9± 


and (64) for the spin components in (47) era 

(65a) 
(656) 
(65 c) 

It follows immediately that g_^_ ■ g_ = 0, 
f9- = if X g^ and f ■ g_ = f ■ g_ = 0. According to ( |6^ , g_ = 0 when / || /. 
Moreover, from (656)-(65c) it follows that g_^_ = g_ = 0 when / = / = 0 . The 
Eilenberger equations read 

Iihvf\/jif + {e-p-)f = u+■ f - gA-g^ ■ A ( 66 a) 

Iih{vfVn)f + {e-i^-)f- f 

= - g+^ - gA- ig_ X A, (666) 

which must be solved together with (65a), (656) and ( |65c ). These equations have to be 
complemented by self-consistency equations for the self-energies. Near Tc all off-diagonal 


quantities (/, /, A, A) are small, and then from (65a)-(65c) follows that g^ can be 


neglected being second order in /, /, and g can be replaced by its normal state value 
[—i 7 rsign(£„) in Matsubara representation]. If one introduces the matrices 


(67) 




^+,x 

^+,y 


\ 

I/= 

^+,x 

Z/_ 

W-,z 

r 

1 • 



^+,y 






^+,Z 

W_^y - 

-iT^-,x 

V- 


and 







^ 9 

9+,x 

9+,y 

9+.Z 

\ 

G= 

9+,x 

9 

i9-,z 

~'^9-,y 




-ig-,z 

9 

^9—,x 



\ 

i9-,y 

1 

1 

9 

/ 

as well as the vectors 






( f ] 


(^ \ 



F = 

fx 

, A = 

^x 




fy 


Ay 

"s 



^ fz ) 






( 68 ) 


(69) 
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then the system 


-( 666 ) can be written compactly as 


2 £‘ — 2 ^)F — ^ G 


(70) 


The four eigenvalues of 2{e— V)/fiVf determine the wavevectors in the system, and can 
be calculated as 






+ Pi 


/(hvf). 


(71) 


2{e — P-) ± 

For a singlet superconductor in proximity contact with a ferromagnet with exchange 
field J, assuming for the following discussion p_ = 0, = 0, = — J, A = 0, the 

Eilenberger equations simplify to 

(lihvfVji + e) f = - J ■ f - gA (72a) 

(^^ihvfVn+e) f = -Jf-g^A. (726) 

where we have introduced e = e + iag for Im(£) > 0 and e = e — ia^ for Im(£) < 0 
to account for possible spin-flip scattering with rate cks = h/rg, and Tg is the spin-flip 
scattering time. The wavevectors are then given by 

^ 1,2 = 2(e ± \ J\)/hvf, 6 ) 3^4 = 2e/hvf. (73) 

In a superconductor-ferromagnet hybrid structure, the superconducting order parameter 
vanishes in the ferromagnetic regions, A^f = 0. The proximity effect manifests itself 
in nonzero pair amplitudes /, / 7 ^ 0 in the ferromagnet. In the superconductor, the 
exchange field vanishes, J = 0, expressing the non-coexistence of superconducting and 


ferromagnetic orders. Equation (73) defines the propagation and decay of Cooper pairs 
into the ferromagnet. If one replaces e by the first Matsubara energy ink^T, then one 
obtains an exponential decay on the length scale = hvf/2(TrkBT+as). In addition, two 
of the components oscillate on the length scale = hvf/2\J\ (with wavelength 27r^j) 
along the direction Vf. After averaging over all directions, this leads to an algebraic 
decay ~ l/x in ^-direction away from the interface, before on a larger length scale 
exponential decay sets in HSU. The components with wavevectors ki^ 2 , oscillating on 
the scale ^j, correspond to eigenvectors with / = ±/J/|J| (i.e. f X J = 0), and 
represent an equal weight superposition of the singlet amplitude (S = 0,3^ = 0) and 
the triplet amplitude with zero spin-projection to J (S = l,Sz = 0). The (degenerate) 
components with wavevectors 6 : 3 ^ 4 , insensitive to J, correspond to an eigenvector where 
/ = 0, / • J = 0, and represent equal-spin pairing amplitudes with respect to the 
quantization axis J (S = 1, Sz = ±1). For increasing exchange field, the pre-factor 
for the proximity amplitudes oscillating with wavelength is of order A/| J|, and thus 
vanishes in the limit of large exchange splitting. 


4 . 5 . Usadel equations 


In the dirty (diffusive) limit, the transport equation (48) can be greatly simplified. A 


strong scattering by impurities averages the quasiclassical propagator over momentum 
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directions. The Green’s function may be expanded in the small parameter tTc (r is the 
momentum relaxation time) following the standard procedure |218[ 1186) 

g{Pf,R,en) g^^\R,£n) + g^^HPf,R,£n) (74) 

where the magnitude of g^^'> is small compared to that of g^^\ The impurity self-energy is 
related to an (in general anisotropic) lifetime function T{pf,Pf) |186j . Substituting (74) 
into ( [ 4 ^ and (49), multiplying with Nf{pf)vfj{pf)T{p'j:,Pf), averaging over momentum 
directions, and considering that Jr/h <C 1, one obtains (with the help of the 
normalization condition) {Nf{pj)vfj{pj)g^^\pj:))p^ = NfJ2k{Djk/iT^)g^^'''^k9^^\ where 


D, 


jk 


N]{ 




Pf' Pf 


(75) 


is the diffusion constant tensor. The function g^^'* obeys the Usadel transport equation 
(hereafter we omit the superscript (0) and apply the Einstein summation convention) 

V, [gVfg) = 0 (76) 


£T3 - ho, g 


+ 


TT 


where ho = (-h^/(P/)[7ext(P/) + hmf(P/)])pj/-h^/, together with the normalization condition 

(77) 


f = 




The diffusion constant is a material property. It is spin independent in agreement 
with the approximation of treating all exchange energy effects as perturbation. In the 
diffusive limit, the different components of g and / in the spin space are related through 
(65a), (656) and (65c), since the Green’s function g averaged over momentum fulfills the 
same normalization condition as before averaging. A vector potential enters in a gauge 
invariant manner by replacing the spatial derivative operators by (see e.g. m) 


v,a: ^ ax = - i 


ThAi,X 

h 


The charge current density is obtained from g{R;£n) via 

2eNfDkj 


3k — 


TT 


[/•£»,/-/• a,/], 


Sn 


with dj = Vj — and the spin-vector current density in spatial /c-direction 

2eNfDkj 


Jk — 


TT 


V,g + r xa,/]. 


(78) 


(79) 


(80) 


Note that from Eqs. (656)-(65c) follows that g = [Im(/*/) + iRe{f) x Im(/)]/Im(5'o) 
at Matsubara frequencies. As for the case of Eilenberger equations, an important part 
of the problem is the formulation of appropriate boundary conditions (2201 ion]. In 
quasiclassical theory this is usually a highly non-trivial task, as boundary conditions 
are in general non-linear. A powerful way to implement boundary conditions and to 
solve Usadel equations is also here the Riccati method [721 1222) 1223) . Modern versions 
of boundary conditions for Usadel equations are found in Ref. nsa for weak spin 
polarization, small spin-dependent scattering phase shifts, and small transmission. A 
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general formulation appropriate for arbitrary transmission, spin polarization, and spin- 
dependent phase shifts has been derived in Ref. |224) . The tunneling limit of this 
boundary condition has been used in Ref. HM], and was independently introduced in 
Ref. [Hsn]0 

Near the critical temperature Tc, the pair amplitudes / aud f are small aud the 
Greeu’s fuuctiou g deviates ouly slightly from its value iu the uormal state, so that the 
Usadel equatious cau be liuearized aud take the simpler form [731 [2261 [ZH 


f — — J • / + *vrA (81) 

+ / = - J/, (82) 

where e is defined as after Eq. ( 726| ). The characteristic wavevectors for isotropic systems 
{Dij —>■ DSij) are now given by (we chose the roots such that Im(/c) > 0) 

kl,2 = (1 + 0 '\/(^ i ^3,4 = (1 + . (83) 


If we again substitute e —>■ this leads to exponential decays for all four solutions. 

However, there are two long-range solutions with decay length 

and two short-range solutions showing damped spatial oscillations with decay length 


y = I __ 

' \ y(7rfcBr + a,y + I jp + (tt^bT + a,) 

aud iuverse oscillatiou wavevector 

f ^ _ 

' V y(7rfcBr + a,y + I J|2 - (tt^bT + a^) 

Note that < ^2 except when T = 0 and H/ts = 0, thus the oscillations are strongly 
damped, so that one cannot expect to observe many oscillation periods. This is in 
contrast to the clean case, where the oscillation wavelength is entirely decoupled from the 
decay length, which can become very long for low temperatures. Thus, many oscillation 
periods are expected to be observable in experiment. For the diffusive case, on the other 
hand, both length scales are temperature dependent. This allows for oscillatory behavior 
as a function of temperature, an effect absent in the clean limit. For | ^ k^T, ag both 

length scales merge and approach [227] 

With iucreasiug exchauge splittiug, the decay leugth shriuks to zero aud the short-rauge 
proximity amplitudes are expected to be suppressed to uumeasurable size, uuless oue 
fiuds a mechauism to create uew types of proximity amplitudes. 

A recently suggested alternative boundary condition |225] assumes the simplified case of coinciding 
spin-quantization axes for reflection and transmission matrix elements of each scattering channel within 
a channel-conserving approximation, and has a considerably narrower range of applicability. 
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4-6. General features of the S/F proximity effect 


It can be noticed from the transport equation (726) or (82) that the triplet vector / obeys 


in the ferromagnet an inhomogeneous differential equation, implying that / is necessarily 
non-zero if the singlet component / penetrates in the ferromagnet. For spatially constant 
exchange field J, the triplet vector aligns with J. The singlet component and the triplet 
component with the zero spin-projection on J coexist always in the ferromagnet near 
the S/F interface. As they both involve electrons from both spin bands, both are 
characterized by short-range penetration lengths in the ferromagnet. 

On the contrary, if the triplet vector / is non-collinear with J, triplet components 
with nonzero spin-projection on J are produced. Since these correspond to equal- 
spin pairing, they are not limited locally by the paramagnetic interaction with the 
local exchange field and have long-range scales in the ferromagnet. A misalignment 
between the triplet vector / and the moment J occurs in presence of sudden changes 
in orientation of J. The reason is that / obeys a differential equation and its variations 
in orientation have thus to be relatively smooth. 

The counterpart for the production of triplet components (/ f 0 ) in the pair 
amplitudes is the presence of gf 7 ^ 0 in the diagonal components of the Green function 
Em. As a direct consequence, the density of states for the up and down spin projections 

This feature of the S/F proximity effect has been found 
In the presence of long-range triplet components. 


differs [see Equation (63) 
numerically as early as in 1999 
the particle-hole diagonal Green function components contain also off-diagonal terms in 
the spin-space, a signature of a spin-flip scattering process. 

As a result of the spin splitting in the density of states generated by the S/F 
proximity effect, a spin magnetization SM is also induced near the S/F interface. This 
magnetization leakage has been investigated in Ref. |229) within a model considering a 


fixed exchange field. In the dirty limit the spin magnetization induced by the proximity 
effect is given by |186[ 1151] 


SM{K) = 2NfksT Re [g (il, £„)], 


(87) 


with hes— (i—Since the triplet vector / is also induced in the superconductor 
near the S/F interface via an inverse proximity effect, the vector g characterizing the 
magnetic correlations penetrates also in the snperconductor according to the relation 


(65c), which for the diffusive limit reads Im(go)Re(^) =Im(/*/) (all functions of £^). 

It is possible to convince oneself that the sum over Matsubara frequencies in the 


expression (87) is nonzero. Indeed, as noticed in Refs. [nnEsoiEsi], in the diffusive 


limit the triplet components are odd functions of the Matsubara frequencies while 
the singlet amplitude / is an even function of £„. Note that in the diffusive limit 


9o(-£n) = -9o(£n), s(-£n) = g(^n) ■ 


( 88 ) 


For unitary pairing states, / x / = 0 , the relations (|656[)-(|65c[) between the vectors 


g and / get simplified. In this case, it can be seen that if the gap A can be chosen real 
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(and thus the singlet amplitude f{en) is real), then the triplet vector f{en) is purely 
imaginary, i.e. f = f* = —f (taking into account that go = —go, i.e. go is purely 
imaginary). As a result, one obtains the simpler relations 


g = 9, 9o9 = // 

These simplihcations arise also when / 


J 1221]. Combining the relations 


(89) 

and 


(88), it follows that the spin-vector part g is an even function of in the diffusive limit, 
which demonstrates that the induced spin magnetization SM is in general non-zero. 

Near T^, the singlet amplitude / and the triplet vector amplitude are small so that 
g and thus 5M appear to be second order terms [see Equation (65c)]. Accordingly, the 
induced spin magnetization SM penetrating the superconductor, which is negligibly 
small near Tc, increases significantly by reaching temperatures well below Tg. 


4 . 7 . Strongly spin-polarized ferromagnets 

In the case that the spin splitting of the spin bands in a ferromagnet is much larger 
than the superconducting gap in the adjacent superconductor, the above mentioned 
theories must be modified in various respects. First, Usadel approximation assumes 
that the inverse life time of quasiparticles due to impurity scattering averaged over 
the Fermi surface, h/r, is much larger than all energy scales appearing in the Usadel 
equation, including the magnitude of the exchange field \J\. If this is not the case, the 
spatial variation of superconducting correlations into the ferromagnet happens on length 
scales comparable to or shorter than the mean free path, and Usadel theory cannot be 
used within the ferromagnet on the same footing as the exchange field. One must 
resort to Eilenberger equations in this case. If the exchange splitting is, however, not 
much smaller than the Fermi energy scale, then even Eilenberger-Larkin-Ovchinnikov 
theory cannot accommodate the exchange field in the way described in the previous 
sections. A modified version of Eilenberger and Usadel equations can be used for the 
case that \J\ ~ Ep- In this case, there exist two well separated fully spin-polarized 
Fermi surfaces in the system, and the length scale associated with h/lpj^ ~ Pfl\ ^luch 
shorter than the coherence length scale in the ferromagnet. Superconducting pairs of 
the type tt 4-4- c'-re still long-ranged in such a system, however pairs of the type 
44 and 4-t c'-re negligible within quasiclassical approximation. Both Fermi velocity and 
density of states at the Fermi level are spin-dependent. The same holds for diffusion 
constant, and coherence length. The quasiclassical propagator is then spin-scalar for 
each quasiclassical trajectory in the ferromagnet, 

and only has a particle-hole discrete degree of freedom. Similarly, all mean field self 
energies have the same structure (spin-scalar, only 2x2 particle-hole matrices), and the 
Fermi liquid interactions are replaced by spin-scalar interactions, which in the simplest 

■fA ^ 

case do not mix the spin bands: A®, A“—>■ Ay^, U®, Vy, (in a more general 
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case, the Fermi liquid interactions could be of the form Eilenberger 

equation and Usadel equation have the same form as before for each separate spin band. 
The spin-resolved current densities are given in the ballistic case by 

j't = h = , (91) 

and the spin magnetization by 

Mz = Mn{Bz) + /cbE ■ (92) 

The expressions in the diffusive limit are modihed accordingly: e.g. for the spin-resolved 
current density one obtains 

, (93) 

and analogously for spin down. 

5. Pair amplitude oscillations and 7 r-Josephson junctions 

5.1. 0 — TT transitions in Josephson junctions 

Proximity-induced pairs in a superconductor-ferromagnet bilayer are subject to the 
FFLO effect, leading to spatially oscillating pair amplitudes in the ferromagnetic regions. 
These act back on the superconducting singlet pair potential in the superconducting 
region. How exactly this interaction between pair potential and FFLO amplitudes 
takes place depends on details of the interfaces and must be determined numerically, 
however for a sufficiently thin superconducting layer the oscillating nature of the FFLO 
amplitudes ultimately leads to oscillations in the modulus of the pair potential of the 
superconductor (and thus in the transition temperature) as function of the geometric 
dimensions of the hybrid structure. In addition, the oscillatory nature of the FFLO pairs 
in the ferromagnet affects the properties of superconductor-ferromagnet-superconductor 
Josephson devices, leading to the possibility of oscillations between junctions with a 
phase difference of zero and junctions with a phase difference of tt in the ground state 
when a certain control parameter is changed. 

The idea of Josephson junctions with a tt phase difference in its ground state when 
a magnetic impurity is inserted in the junction was introduced by Bulaevskii, Kuzii, and 
Sobyanin in 1977 |232) . Similar ideas had also been proposed by Kulik in 1965 |233) . In 
1982 it was also shown in a classical paper by Buzdin, Bulaevskii, and Panyukov that 
the pair amplitude oscillates in ferromagnets in contact with a superconductor |234) . 
Further early studies of this effect followed in Refs. |235[ 123611237] . One technological 
problem that became evident quickly was that weakly spin-polarized systems, like 
ferromagnetic Cu-Ni or Pd-Ni alloys, are better suited to observe 0 — tt oscillations, 
as otherwise the proximity amplitudes become so short ranged that extremely thin 
layers are necessary. Typical ferromagnetic alloys that have been used are Cui_a;Nia; 
(CuNi) and Pdi-^iNL (PdNi). The real break-through came in the beginning of the 
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Figure 13. (a) Temperature dependence of the critical Josephson current in a 

Nb/Cu 52 Ni 48 /Nb trilayer device, shown in the inset, for a ferromagnet layer thickness 
d = IS nm. (The lines correspond to theoretical models with constant (dotted 
line) or temperature dependent (solid line) exchange energy.) From Sellier et al 
[243] . Copyright (2003) by the American Physical Society, (b) The ferromagnet-layer 
thickness dependence of the critical current density for Nb/Cuo. 47 Nio. 53 /Nb junctions 
at temperature 4.2 K. Open circles represent experimental results; solid and dashed 
lines show model calculations. The inset shows a schematic cross section of the SFS 
junctions. From Oboznov et al. [223. Copyright (2006) by the American Physical 
Society. 


2000’s with the experimental verification of the switching between 0 and vr-Josephson 
junctions based on the effect predicted theoretically before. Pioneering work was done 
by Ryazanov and co-workers who studied transitions between zero- and yr-states as 
function of temperature [2381 I239[ I240[ I227j and thickness of the ferromagnet layer 
[24011227| in Nb/CuNi/Nb trilayers (see Figure [I^. Closely following were experiments 
by Kontos and co-workers who studied the zero-7r transition as function of barrier 
thickness in Nb/Al/AROs/PdNi/Nb junctions |24f ) . by Blum et al. [242] . who use 
Nb/Cu/Ni/Cu/Nb junctions and vary both temperature and thickness of the Ni layer 
(thus employing a much stronger spin-polarized ferromagnet than the ferromagnetic 
alloys that were used by the other groups), and by Sellier et al. studying a temperature 
induced 0 — tt transition in Nb/CuNi/Nb (see Figure [I^, as well as the appearance of 
half-integer Shapiro steps [24311244] . 

The critical Josephson current density as function of layer thickness in the limit 
fcsT, tts \J\ is described in the clean limit by formula mi 





(94o) 


where do is a fit parameter which is usually assigned to a “dead layer”,i-e. a layer 
of suppressed magnetism at the interface. In the diffusive limit the corresponding 
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Figure 14. (a) Characteristic voltage of a Nb/Co/Nb Josephson junction as a 

function of Co thickness at 4.2K. The dashed line is a fit to theory. Inset: A focused 
ion beam micrograph of a typical Nb/Co/Nb Josephson junction. From Robinson et 
al. PT7] . Copyright (2006) by the American Physical Society, (b) Critical current 
vs applied magnetic field obtained for Nb/Cu/Co/Ru/Co/Cu/Nb circular Josephson 
junctions with total Co layers thickness of 6.1 nm. From Khasawneh et al |250| . 
Copyright (2009) by the American Physical Society, (c) jc{dF) dependence for an SFS 
junction with F=Ni. Two fits appropriate for the dirty limit theory for data from 
regimes I+II+III and just II+III are shown as dashed lines, (d) Ic{T) dependence 
for SFS Josephson junctions. A temperature induced 0 to tt transition is observed 
(sample dp = 2.91 nm). (c) and (d) from Bannykh et al. |249| . Copyright (2009) by 
the American Physical Society. 


expression for long (c?f S> Ci) junctions is 

(ip — do 


Icidp) ^ e cos 


6 


(946) 


with the length scales ^ 1^2 as in (85) and (86). A study by Pugach et al. |245) bridging 
the two limiting cases showed that the two length scales and ^2 may exhibit a non¬ 
monotonic dependence on the properties of the ferromagnetic layer such as exchange 
field or electron mean-free path. 

Further studies of strongly spin-polarized interlayers were performed, using various 
materials, among those Ni in Ref. |246) . Co, Ni, and Py (Ni 8 oFe 2 o) in Ref. |247) . and Py 
(Feo. 75 Coo. 25 ) barriers in Ref. [248] . A detailed study, presented in Ref. [249] . with Ni 
interlayers using the formulas (94a), (946) above has been done recently, indicating that 
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these structures are in the dirty limit. The thicknesses of the Ni layers ranged between 


1 and 5 nm. In Figure 14 examples for 0 — vr oscillations in a Nb/Co/Nb structure are 
shown. In order to eliminate unwanted effects due to the variation of the electromagnetic 
vector potential across the junction due to the intrinsic magnetic flux in the junction, a 
special geometry was used where two Co layers were exchange coupled by a thin Ru layer 
in between to align antiferromagnetically. This cancels all the unwanted effects to a high 
degree, and a previously wildly fluctuating critical current vs. magnetic field curve turns 
into an excellent Fraunhofer pattern when the Ru interlayer present, as seen in Figure 
14 (b) |25()) . With this improvement high-quality 0 — tt transitions are observed both 


as function of temperature and of ferromagnet thickness. These results complement the 
results shown in Figure [T^ which are for ferromagnetic alloys. The main difference is 
a one order of magnitude smaller ferromagnet layer thickness necessary to observe the 
effect, which is an experimental challenge. 

The effect of additional normal and/or insulating layers between the superconductor 
and the ferromagnet was studied theoretically in detail in Ref. ISH. It was found 
that even a thin additional normal conducting layer may shift the O-vr transitions to 
larger or smaller values of the thickness of the ferromagnet, depending on its conducting 
properties, and for certain parameter ranges a O-tt transition can even be achieved by 
changing only the normal layer thickness. 

A promising setup are so-called double-proximity structures, in which two 
superconductors are connected by a weak link of a normal-metal/ferromagnet bilayer 
or a ferromagnet/normal-metal/ferromagnet trilayer forming a bridge between the 
superconducting banks. The proximity effect acts here twice: first to provide pair 
amplitudes from the superconducting banks into the bilayer or trilayer, and second 
at the interface between the normal metal and ferromagnet. Such structures have 
been proposed by Karminskaya and Kupriyanov |252) . and experimental realization 


in terms of hybrid planar Al-(Cu/Fe)-Al submicron bridges is reported in Ref. [253]. 
The oscillation periods and damping lengths of the critical Josephson current can 
be much longer in such devices than in the conventional setup. Various geometries 
have been theoretically studied in subsequent work, including (SN)-(NF)-(SN), (SNF)- 
(NF)-(SNF), (SNF)-N-(SNF), and S-(NF)-S junctions, in order to optimize practical 
performance |251|. 


5.2. Phase sensitive measurements 

The current-phase relation (CPR) as a function of temperature in a Nb/CuNi/Nb 


junction was measured by Frolov et al. (see Figure 15) |255j . Using an rf SQUID, 
the change of the current-phase relation from zero to 7r-j unction behavior is clearly 
visible. This allows to observe the sign of the supercurrent. A 7 r-junction is a 
Josephson junction with a negative critical current A, as its current-phase relation 
(CPR) is Ij{(p) = I Jc| sin( 9 ?-|- tt) = —|Id sin( 9 ?). As conventional measurement of the 
current-voltage characteristic of a Josephson junction is not sensitive to the sign of the 









48 



Figure 15. (a)Circuit for measuring the current-phase relations of a SFS junction, 
(b) Current-phase relation derived from rf SQUID modulation curves showing the 
transition to a tt Josephson junction as the temperature is lowered, (a) and (b) from 
Prolov et al. |255j . Copyright (2004) by the American Physical Society, (c) A drawing 
of the cross section of a single SQUID junction, (d) Critical current modulations 
showing the expected no-shift between a 0 — 0 SQUID and a tt-tt SQUID and a shift 
of 4>o/2 between a 0 — tt SQUID and a 0 — 0 SQUID. As the of the 0 — 0, 0 — tt, 
and TT-TT SQUIDs are different the Ic{B) curves are observed at different temperatures 
in order to measure in the same range of critical currents for each couple, (c) and (d) 
from Guichard et al. |257) . Copyright (2003) by the American Physical Society. 


supercurrent, it was necessary to include the Josephson junction in a multiply connected 
geometry. The sign was in Ref. |255) observed in an rf SQUID configuration [see Figure 


a)] by shorting the electrodes of the junction with a superconducting loop. If that 
loop contains a tt junction in zero external field it will exhibit a spontaneous circulating 
current, generating a flux of 'Fo/2 which can be detected by a SQUID magnetometer or a 
Hall probe. In Figure [T^ (b) the transition between a zero and a tt state is unambiguously 
observed (although residual magnetic fields leading to shifts in the CPR curves did not 
allow to pin down if it was a O-tt or a tt-O transition). 

Phase sensitive measurements of the ground state of ferromagnetic Josephson 
junctions using a single dc SQUID have been performed by Ryazanov and co¬ 
workers |256j for Nb/NbOj^/CuNi/Nb junctions, and by Guichard et al. |257j for 


Nb/NbOjy/PdNi/Nb junctions, showing that the sign change of the Josephson coupling 
is observed as a shift of half of a flux quantum d>o = ^ in the SQUID diffraction pattern 
(see Figure 15). A 0-0 and a tt-tt junction show identical critical current modulations 
with flux, whereas a O-tt junction exhibits a pattern shifted by $o/2 with respect to 
that for 0-0 and tt-tt junctions. In this experiment a single dc SQUID with two SFS 
junctions embedded. The ferromagnetic layer thickness dpi and dF 2 can then be chosen 
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cross-section 
of junction 


Figure 16. Temperature dependent magnetic flux produced by a tt- and a 0-loop 
when cooling (a) in zero field and (b) in a magnetic field equal to half a flux quantum 
in the loop. Below T ^5 K the loop develops a spontaneous current when cooling 
down in zero field, while the magnetic flux through the 0—loop remains zero. When 
applying a field equal to half a flux quantum the roles of the tt- and the 0-loops are 
interchanged. At low temperatures the spontaneous flux saturates at a value close to 
^o/2 in each case, (c) In zero external magnetic field, a spontaneous flux ^ = ^o/2 
is required to maintain the condition of fluxoid quantization if a 7r-j unction is inserted 
into a superconducting loop (d) Cross section of the planar Josephson junction. From 
Bauer et ai |258| . Copyright (2004) by the American Physical Society. 


to correspond to zero or tt coupling independently, leading to the possibility of a 0 — 0, 
a TT-TT, a O-TT, or a tt-O SQUID. For a dc SQUID with negligible loop inductance and 
equal junction critical currents the modulation of the critical current with applied flux 
is given by P6] 


^c(®ext) — 2/o 


COS TT 




ext 


$0 


(^12 


(95) 


where ^12 is the sum of the internal phases (0 or tt) in the two junctions of the SQUID. 
Thus, for a 0 — 0 and a tt-tt SQUID the /c($ext) patterns are identical, whereas they 
are shifted by half a flux quantum if 5 i2 = tt, i.e. a 0 — tt or a tt-O SQUID. This shift is 
shown in Figure [l5| (d). 

The presence of spontaneous magnetic moments in snpercondncting (Nb) loops 
containing a ferromagnetic (PdNi) tt junction was experimentally demonstrated by 
Bauer et al. (see Figure [l6|) |258) . The loops were prepared on top of a micro-Hall sensor. 
The authors observed asymmetric switching of the loop between different magnetization 
states when reversing the sweep direction of the magnetic field. The presence of a 
spontaneons current near zero applied field was studied as function of temperatnre, and 
the magnetic moment approached half a flux quantum at low temperatures. 

An rf SQUID geometry where the macroscopic gronnd state of a ferromagnetic 
(PdNi) Josephson junction shorted by a 0 weak link was experimentally investigated by 
Della Rocca et al. and showed spontaneous half quantum vortices with random sign 
(±'I>o/2) [259] : it was fonnd that 0 — tt junctions behave as classical spins. 
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Figure 17. (a) Sketch (cross section) of the investigated SIFS O-tt Josephson junctions 
(Nb/Al-Al 203 /Nio. 6 Cuo. 4 /Nb) with a step in the ferromagnet layer thickness, (b) 
Domains of existence of 0-, tt- and <y9-Josephson junctions. The star shows the position 
of the investigated Josephson junction in (c). (c) Experimentally measured Ic{B) 

(black symbols) and theoretical predictions (red symbols), (a)-(c) from Sickinger et al. 
|2fi8| . Copyright (2012) by the American Physical Society, (d) Ic{B) of 0—junction 
(red triangles), 7r-junction (blue triangles), and O-tt junction (black spheres). Prom 
Weides et al. |2h4| . With kind permission from Springer Science and Business Media. 


5.3. ip-Josephson junctions 

An interesting separate development concerns a geometry where a step in the thickness of 
the ferromagnet is present in a superconductor-ferromagnet bilayer, coupled to another 
superconductor via an insulating barrier. On one side of the step one has a 0 junction, 
on the other side a tt junction. A half flux quantum (semifluxon) is trapped in such 
a structure at the step |260[ 12611 [262], and a supercurrent circulates in the structure, 
similarly as in a 0 — tt SQUID |26,H[ 1264112651 . Such structures were studied in detail 

by Weides, Kohlstedt, Koelle, Kleiner, Goldobin and co-workers. In particular, there 
exists in such structures the possibility of a (/^-junction, which has neither 0 nor tt phase 
difference as a ground state. A method to realize a (p Josephson junction by combining 
alternating 0 and tt parts with intrinsically non-sinusoidal current-phase relation was 
suggested in Ref. [Ml- The realization of an SIFS (/^-junction has been experimentally 
achieved in 2012 (see Figure |268j . Properties of zero-7r junctions which act as a 
Josephson junction with an equilibrium phase difference 0 < ip < tt are discussed in 
Ref. [269] . where the current-phase relation is calculated numerically and in certain 
limiting cases analytically. 

In Refs. [2701 I271j . a planar Josephson junction with a ferromagnetic weak link 
located on top of a thin normal metal film is considered. It is shown that this Josephson 
junction is a promising candidate for the realization of a (y?-junction. 
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6. Long-range triplet supercurrents 

6.1. Length scales for superconducting correlations in ferromagnets 

As detailed in the sections [4^ and [47^ proximity induced superconducting correlations in 
ferromagnets with sufficiently strong exchange splitting | J| come in pairs of short-range 
and long-range amplitudes. The former ones are the singlet amplitude and the triplet 
pair amplitude with zero spin projection, Sz = 0, on the magnetization axis of the 
ferromagnet. These short-range amplitudes decay in ballistic structures algebraically 
with a reduced magnitude of order A/|J|, and in diffusive structures exponentially 
on the length scale = y^hDj\J\. The other two components, called long-range 
amplitudes, decay on the thermal coherence length scale given for ballistic structures 
by hvf/{27rkBT + 2as) and for diffusive structures by ^JhD/{2^^kBT + 2as). They are 
characterized by triplet amplitudes with a spin projection = ±1 on the magnetization 
direction of the ferromagnet. Consequently, they are called equal-spin triplet pair 
correlations. Equal spin pairs do not suffer from having to populate spin-split pairs 
of Fermi surfaces. Both electrons of the pair are situated on the same spin component 
of the Fermi surface. For this reason, they behave like in a normal metal, with the 
coherence length determined by the Fermi surface properties of one spin component 
only (which in general differ for Sz — -\-l and Sz ~ ~1)- The quest for such long- 
range pair amplitudes has a long history, and is one of the success stories of interaction 
between experiment and theory. 

6.2. Experimental “pre-history” 

As in the 1990’s the field of spintronics, i.e. functional nanometer-size devices based on 
spin-dependent transport phenomena, developed rapidly [272], a strong motivation to 
search for a long-range proximity effect in ferromagnets was established. Such long-range 
amplitudes would lead to long-range supercurrents in Josephson devices, and would 
ultimately lead to valuable applications. The ultimate goal is to obtain completely 
spin-polarized supercurrents, which would necessarily have to be triplet supercurrents. 

It was for this reason that in the second half of the 1990’s experimental efforts 
intensified to study proximity effects in strongly spin-polarized ferromagnets. Petrashov, 
Antonov, Maksimov, and Shaikhaidarov found in 1994 that the effect of superconducting 
islands, deposited on the surface of ferromagnetic Ni, can be still seen in the resistance 
of the Ni layer distances greater than 2 /rm (exceeding by more than 30 times what the 
length hvfl\J\ would suggest) away |273j . 

Lawrence and Giordano studied in 1996 the resistance as a function of temperature 
and magnetic field in structures containing a narrow (~ 2;um) ferromagnetic strip 
connecting two superconducting films (In/Ni/In and Pb/Ni/Pb) |274] . They observed a 
magnetoresistance dip as function of magnetic field much too large to be accounted for 
by weak localization effects. In addition they found implausibly long phase coherence 
lengths, inconsistent with the usual superconducting proximity effect. In a later study 
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of Sn/Ni/Sn structures, where the Ni was a narrow wire 40 nm wide), they measured 
an unexpectedly long proximity length of about 50 nm (theory predicted 4 nm for these 
structures) when the Sn/Ni interfaces are clean, however with an oxide layer at the 
interfaces they observed an unexplained re-entrant behavior They also noted that 

their proximity length was of the order of the thermal length y^hD/k^T which would 
appear in a normal metal. 

Giroud and co-workers studied in 1998 Co/Al structures with a ferromagnetic 
Co wire in contact with superconducting A1 They observed below the 

superconducting transition that the Co resistance exhibited a significant dependence on 
temperature and voltage, and the differential resistance showed that the decay length 
for the proximity effect was much larger than expected from the exchange field of Co. 

Petrashov and co-workers found in 1999 a giant mutual proximity effect in 
Ni/Al structures with an proximity-induced conductance on the Ni side two orders 
of magnitude larger than predicted by theory 1277]. Aumedato and Chandrasekhar, 
by performing multi-probe measurements on Ni/Al structures, re-evaluated these 
studies and came to the conclusion that superconducting correlations cannot 
extend into a ferromagnet over distances larger than the exchange length, and 
associated the large changes in resistance seen in the previous experiments with the 
superconductor/ferromagnet interface or the superconductor, which had been measured 
in series or in parallel with the ferromagnet [27H]. The controversy lead to a vivid 
discussion about the existence of such long-range proximity components. 

6.3. Theory of long-range proximity amplitudes 

6.3.1. Spiral magnetic inhomogeneities and domain walls The renewed interest created 
by the lack of understanding of long-range superconducting proximity effects in 
ferromagnets lead in the beginning of the 2000’s to new theoretical developments. 
Pivotal was a series of papers by Bergeret, Volkov, Efetov, and coworkers, who studied 
the effect of a spiral Bloch-type inhomogeneity close to an S/F interface, where 
the magnetization vector rotates along the junction direction. They found that an 
eqnal-spin triplet pair amplitude shows a long-range penetration into the ferromagnet 
|ldl[ 1721 1222L127911280] . The authors also noted that this, as an isotropic (“s-wave”) 
triplet pair component, must be a realization of odd-frequency pairing. As impurities 
suppress all anisotropic pairing components, this odd-freqnency amplitnde is the only 
superconducting pairing amplitude present in the ferromagnet. Shortly after that, 
Kadigrobov, Shekhter and Jonson found a similar effect |281j . A realization of these ideas 
motivated by experiment |282j was studied in Ref. |28dj , where the focus of study was a 
setup of a Ho/Co/Ho trilayer sandwiched between conventional s-wave supercondncting 
leads. Holminm is a conical ferromagnet with an intrinsic spiral strnctnre. In Ref. |284j 
a bilayer consisting of an ordinary singlet superconductor and a magnet with a spiral 
magnetic strncture of the holmium or erbium type was investigated by solving self- 
consistently Bogoliubov-de Gennes equations. A re-entrance behavior as function of 
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temperature is observed for certain parameter ranges. 

Another type of magnetic structure discussed in the literature is a domain 
structure of the Neel type, in which the magnetization vector lies parallel to the 
interface and rotates along a direction parallel to the interface. A setup with a chiral 
ferromagnet, exhibiting a homogeneous cycloidal spiral, placed between two conventional 
superconductors was studied in Ref. |285) . Depending on the spiral wave length, O-tt 
transitions can be induced. In the case of a uniformly rotating spiral, however, it 
was shown that only short-range components exist |285] . In a more general case, with 
magnetic domains separated by Neel walls, long-range triplet components are present 
and arise at the domain walls, decaying inside the domains [28611287] . 

Refs. |288l 1289] investigate diffusive SF and SFS structures with a domain wall 
in the F layer, using a Riccati representation of the non-linear Usadel equations ra 
and including spin-mixing parameters for interfaces. They study local density of states, 
induced magnetization, and spin-polarized Josephson currents. It is found that the spin 
polarization of the spin current in SFS structures is determined by the magnetization 
profile and that the spin current shows discontinuities at the zero-yr transition points 
of the critical Josephson current. Similar as in Ref. m, a non-zero spin-current even 
for zero superconducting phase difference is reported. In Ref. fm\ the contribution of 
domain walls to the Josephson current through a ferromagnetic metal is examined for 
both ballistic and diffusive systems. It is found that in the clean limit domain walls 
enhance the Josephson current even for a collinear magnetic domain structure, whereas 
in the diffusive limit a non-collinear domain structure is necessary to enhance the effect. 
In Ref. |291j the influence of the location of a domain wall within a Josephson junction on 
the ground state properties of the junction was the topic of investigation. The authors 
considered both the diffusive and ballistic limits. They find that the location of the 
domain wall determines if the junction is in a zero-state or in a vr-state and influences 
the transition temperature of the junction. 

A different type of system with inhomogeneous magnetic structure was studied in 
Refs. |292[ 1293] . where a ferromagnetic vortex in a mesoscopic ferromagnetic disc is 
brought in contact with a singlet superconductor. Proximity-induced long-range triplet 
components are generated by the ferromagnetic vortex under certain circumstances, 
leading to zero- or yr-j unction behavior depending on the contact geometry. 

6.3.2. Multi-layer geometries Instead of considering inhomogeneous magnetization 
profiles within a ferromagnetic layer as sources of long-range triplet correlations, one 
can also obtain the same effect in a multilayer geometry with non-collinear arrangement 
of the magnetizations of the various layers. In Ref. |294j a diffusive SFIFS junction 
is investigated by self-consistently solving Usadel equations with appropriate boundary 
conditions. It is found that for antiparallel ferromagnet magnetizations the critical 
current is enhanced by the exchange energy, and for parallel magnetization the junction 
exhibits a transition to a yr-state. A switching behavior between zero- and yr-state can be 
observed when going from the antiparallel to the parallel alignment. In Refs. 
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a multilayered superconductor-ferromagnet structure with non-collinear alignment of 
the magnetizations of the ferromagnetic layers is considered, and the resulting long- 
range odd-frequency triplet condensate amplitude is calculated. It is found that this 
setup allows for a Josephson effect with possible zero- or vr-junction behavior. It 
is also pointed out that the chirality of the magnetization profile matters for the 
direction of the Josephson current. In Ref. |295) a fully self-consistent calculation 
was performed for a ferromagnet-superconductor-ferromagnet trilayer setup with non- 
collinear magnetization, and the induced spin magnetization in the entire structure 
was determined. This work includes self-consistency of the order parameter in the 
superconducting layer at arbitrary temperatures, arbitrary interface transparency, and 
any relative orientation of the exchange fields in the two ferromagnets. The long-range 
Josephson effect through a ferromagnetic trilayer was the topic of Ref. |296] . where 
conditions were derived under which the long-range triplet proximity effect dominates 
the short-range proximity effect. Both diffusive and clean limits were studied. 

Halterman and Vails have investigated numerically ballistic S/F/S trilayer and 
S/F/S/F/S pentalayer Josephson junctions by self-consistently solving Bogoliubov- 
de Gennes equations, concentrating on the thermodynamic stability of zero- and tt- 
j unctions as function of material parameters and geometry [2nzi[2nH]. 

Zero-TT transitions in Josephson junctions with a diffusive pentalayer-structure of the 
form F/S/F/S/F with misaligned magnetizations in the three ferromagnetic layers were 
explored in Ref. |M?| . The transition temperatures for zero- and yr-j unctions as function 
of misalignment angle were determined numerically using an effective and fast procedure, 
which is described in detail in this publication. In Ref. PnOl a diffusive S/F'/FF/F'/S 
structure is discussed, with two different ferromagnetic materials F and F'. The middle 
FF double layer is considered to be parallel or antiparallel magnetized, whereas the 
F' layers were allowed to be non-collinear. This setup is in particular adapted to the 
experiments by Khaire et al. IM]. A similar study for clean or moderately diffusive 
materials was performed in Ref. |3()2) . 

An SF 1 F 2 S junction in the ballistic case (including moderate disorder in the 
ferromagnets) was considered in Refs. |d(),H[ 1504] . with misaligned ferromagnetic layers 
Fi and F 2 . It was found that the long-range spin-triplet correlations lead to a dominant 
second harmonic in the Josephson current-phase relation of asymmetric junctions 
an effect also found for diffusive structures in Ref. IHOi. In the latter reference this 
phenomenon is traced to the long-range coherent propagation of two triplet pairs of 
electrons, a process first pointed out and discussed in Ref. 11721. and called there “crossed 
pair transmission”. 

Magnetic moment manipulation by triplet Josephson currents in Josephson 
junctions with multilayered ferromagnetic weak finks is discussed in Ref. [HOT]. It 
is reported that by tuning the Josephson current, one may control a long-range 
induced magnetic moment, which appears due to non-collinear magnetization of the 
layers. Alternatively, applying a voltage one can in such a geometry generate an 
oscillatory magnetic moment. The spin-switching behavior as function of rotation of 
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the magnetization in one ferromagnetic layer in diffusive SFSFS and SFSFFS Josephson 
junctions was investigated in Ref. |3()8) . 

The proximity effect in clean superconductor-ferromagnet structures caused by 
either the spatial or momentum dependence of the exchange field was discussed in 
Ref. |309j . It was shown that in symmetric junctions the long-ranged proximity effect 
is present for any non-collinear ferromagnetic moment orientation and a dominant 
second harmonic appears for any asymmetric junction. The Josephson coupling between 
two s-wave superconductors separated by a ferromagnetic trilayer with non-collinear 
magnetization was reconsidered in Ref. EIDl. where in particular the dependence on 
strength of the exchange field and thickness of the interface layers was studied, and the 
competition between long-range and short-range components was investigated. 

The appearance of a spin supercurrent in a diffusive SF 1 F 2 S Josephson junction was 
subject of study in Ref. iHm. It was found that a spin current with spin polarization 
normal to the magnetization vectors flows between the ferromagnetic domains. This spin 
current is even as function of superconducting phase difference, and odd as function of 
misalignment angle between Fi and F 2 . 

A diffusive S-(FNF)-S double proximity structure, in which an FNF trilayer with 
mutually misaligned ferromagnetic magnetization vectors bridges two superconducting 
banks, was investigated in Ref. |312j . In this case long-range equal-spin triplet 
correlations are generated by the misalignment of the two ferromagnetic layers in the 
bridge. 

6.3.3. Singlet-Triplet conversion at interfaces The mechanisms discussed in the 
previous subsections employ a quasiclassical approximation, which is not valid for length 
scales as small as the Fermi wavelength. In fact, for the mechanism to work it is crucial 
that short-range components can enter the ferromagnet over a sufficiently long distance 
to feel the magnetic inhomogeneity, which was assumed to vary on the scale much larger 
than the Fermi wavelength. This poses the problem of what happens at interfaces 
between a superconductor and a ferromagnet with large exchange splitting, so that the 
length scale becomes comparable to the Fermi wave length. 

An appropriate mechanism for structures with strong exchange splittings was 
proposed, and first exemplary studied for the extreme case of a half-metallic ferromagnet 
(where one spin band is insulating and one spin band metallic), in [1551 11641 11071 
EE]- The mechanism is based on spin-dependent scattering phases under reflection 
and transmission, which in combination with a misalignment with respect to the 
ferromagnetic bulk magnetization of the interface magnetic moment in an atomically 
thin interface layer leads to the creation of long-range triplet pairs. It was shown 
subsequently, that this mechanism persists for the full range from ballistic to diffusive 
systems ina. Note that such a mechanism invokes the magnetic inhomogeneity directly 
in an interface harrier, rather than in a metallic region. It thus differs from the 
mechanism suggested by Bergeret et al. Indeed, the two mechanisms are complementary, 
with the first one being important for ferromagnets where the exchange splitting is 
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even frequency components odd frequency components 



x/| x/^ 

superconductor half metal superconductor superconductor half metal superconductor 



Figure 18 . (a) All symmetry components in a ballistic S/F/S Tr-Josephson junction 

with a half-metallic ferromagnet (H) as F-layer (self-consistent calculation for model 
barrier as in Fig. 10 (b), with Vs = 0.5, d = 0.5//co); barrier spin-polarized in x- 
direction, half-metal spin-polarized in ^-direction; temperature T = 0.2Tc; (b)-(e) 
Non-monotonic temperature dependence of critical Josephson current density in an 
S/F/S junction, normalized to its value at zero temperature; (b) for half-metallic 
ferromagnet, dependence on the mean free path in H for a fixed junction length 
L = = dvH where vh is the Fermi velocity in H; (c) as in (b), but showing 

the dependence on junction length L for fixed £h = l-5Cc; (d) for a ferromagnet 
in the ballistic limit, with various spin polarizations P = {pfi — Pf^)/{pfi P Pf^)] 
L = hvF/‘2FTc with vf the Fermi velocity in F. (e) IcRn product and normal-state 
resistance RnA [in units of {e^NsVs)~^^ where Ns and vs are density of states per 
spin at the Fermi level and the Fermi velocity in the superconductor, respectively] as 
function of P for T = 0.5Tc, and various strengths of interface transmission. In all 
plots the interface misalignment angle is a = 7r/2. (b) and (c) after |133) . (d) and (e) 
from ina. Copyright (2009) by the American Physical Society. 


small compared to the Fermi energy, and the second one being important for the 
case of comparable exchange splitting and Fermi energy, both being much larger than 
all superconducting energy scales. The crossover between these two asymptotic cases 
cannot be described within either theory, as no useful expansion parameter is present 
in the crossover region. 

The mechanism of long-range triplet pair creation at interfaces was generalized to 
Josephson structures involving strongly spin-polarized ferromagnets with two itinerant 
spin bands in Ref. m- 


6.3.4- Supercurrents in strongly spin-polarized itinerant ferromagnets and half metals 
We discuss first the main predictions of a long-range Josephson effect in a structure 
involving a long half-metallic ferromagnetic region, as discussed in Refs. [miTUTlim 
As shown in Fig. which is for a fully self consistent calculation using a 
model barrier as in Fig. 10 (b), all possible symmetry components shown in Fig. |^are 
indeed generated at the interface. This simply reflects the fact that both spin rotational 
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symmetry and inversion symmetry are broken by the interface. The calculation assumes 
that the effective magnetic moment of the interfaces is perpendicular to the bulk 
magnetization of the half-metallic ferromagnet. As a consequence, according to Eq. 
(246), long-range amplitudes are generated, linking the two singlet superconductors 
[see lower panels in Fig. 18 (a)]. Thus, we have an example of an indirect Josephson 
coupling lEa, via a long-range triplet component which is generated in the first place 
from the singlet component by microscopic processes in the interface region. Note also, 
that the thermodynamically stable configuration for identical interfaces is a 7 r-j unction 
(for asymmetric junctions the appearence of a cpo-junction is predicted |l()7l 1133) ). for 
which the system tries to minimize the odd-frequency triplet amplitudes of Type D in 
the structure. This is a general principle, as odd-frequency amplitudes increase the free 
energy density, instead of decreasing it as even-frequency singlet amplitudes do. Only if 
the p-wave components (and all higher orbital components) are suppressed by impurity 
scattering, can the odd-frequency triplet amplitude dominate. 

A prominent feature in such a structure is a non-monotonous temperature 
dependence of the critical Josephson current density, with a pronounced maximum at a 
temperature which can be linked to the Thouless energy of the ferromagnetic layer (see 
Fig. (c) and Ref. |b33) ). When the junction becomes effectively long compared with 
the diffusive-limit coherence length, L S> = \/Co^h/ 3, the current is dramatically 
suppressed and the peak is shifted to a lower temperature. This maximum moves also 
to low temperatures in diffusive structures, as is seen in Fig. (b), which shows its 
dependence on the mean free path in the ferromagnet. If one replaces the half-metallic 
ferromagnet by a strongly spin-polarized ferromagnet with two itinerant spin bands, then 


the maximum remains only for sufficiently strong spin polarization, as seen in Fig. 18 


(d). If one projects the Fermi surfaces on the contact plane, then only the “half-metallic” 
momentum directions contribute to the maximum, which are inside the Fermi sea for 
one spin projection and outside for the other m Furthermore, there is an optimal 
spin polarization which maximizes the ARn-product of the Josephson junction, which 
for the particular model interface studied in Ref. dza is around P ^ 0.3. The green 
shadowed regions in Fig.[^(d) for small and large spin polarizations are regions beyond 
the validity of the theory, as other processes which are neglected become relevant. 

Following the early theoretical work on triplet supercurrents in half-metallic 
ferromagnets dssi usi, which treats a fully developed triplet proximity effect self- 
consistently, numerous groups confirmed the existence of triplet supercurrents under 
various conditions in half-metallic ferromagnets. In Ref. ra the work of Refs. dsa 
116411107) was generalized to include impurity disorder bridging between the two limits 
of ballistic and diffusive transport. In this work, the full Eilenberger equations with 
impurity self energy in self-consistent Born approximation was solved. In Ref. Ha a 
fully developed triplet proximity effect, just like in Ref. usa for the clean limit, was 
treated for the diffusive limit. Further studies in the diffusive limit were performed 
in Refs. |3T3ll3M[22i] . The role of magnon creation in the half-metallic ferromagnet 
was studied in BEI, with a resulting temperature dependence of the critical current 
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Figure 19. (a) S/F/S Josephson junction with a strongly spin-polarized ferromagnet 

with magnetization M, having interface barriers with misaligned magnetic moments 
mi and m 2 . Only equal-spin pair amplitudes are present in each of the spin bands of 
the ferromagnet. The azimuthal angles of mi and m 2 with respect to the magnetization 
axis M are denoted and (^^ 2 . (b) Critical current in positive (/c+)/negative (/c-) 
bias direction vs A(/9 = ^2 — (c) The zero current equilibrium phase difference 

^Xeq = (X2 — Xi)eq VS A(/9 varies from tt to 0. Here, T = 0.2Tc, P = 0.21. The polar 
angle of mi and m 2 is a = 7r/2. The length of the junction is L = hvs/^nTc^ with vs 
the Fermi velocity of the superconductor, (b) and (c) from |172| . Copyright (2009) by 
the American Physical Society. 


reproducing the temperature dependence found in dna [inzi. The case when the 
superconductor is unconventional was studied in |316l 131711318) . In Refs. |319[ 1215) 
the problem was advanced on the analytical level assuming a constant singlet order 
parameter in the superconductor. Self-consistent solution of Bogoliubov-de Gennes 
equations in the clean limit [320] found results in agreement with Ref. [T551 [TOT] , 
including a strong influence of the junction behavior by subgap Andreev bound states. 

Most of these theoretical studies deal with either the dilfusive limit within the 
Usadel approximation or the clean limit. The full range of impurity scattering from 
clean to diffusive limit for a superconductor-half metal proximity structure was first 
covered in the theory of Ref. BBI. A treatment in the quantum limit was given by Beri 
et al [321] . 

In Ref. [172] a number of effects were pointed out that are absent in the half- 
metallic case, and which are due to the phase coherent transport of Cooper pairs in the 
two itinerant spin bands of a ferromagnet. These are illustrated in Fig. As seen from 


Eq. (246), equal-spin triplet correlations in the ferromagnet acquire an additional phase 
from the configuration of the magnetic moments in the interface region, mi and m 2 , 
with respect to the bulk magnetization M. If the three magnetizations are non-coplanar^ 
there is an important geometric phase involved, which is related to the projection of 
the interface magnetic moments of the two interfaces to the plane perpendicular to the 
bulk ferromagnetic magnetizatioij^ The relative angle between these two projected 
magnetic moments defines the angle A(p = (p 2 — Pi ^ which is a gauge invariant quantity 
independent of the choice of the global spin quantization axis. This angle directly enters 


In the coplanar but non-collinear case there still can be an important geometric phase of tt introduced 
by the triple of magnetic vectors if (mi x M) • (m 2 x M) < 0. 
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the current-phase relation for each spin band, and precisely with opposite sign for the 
two spin projections. Thus, the first-order processes (single Cooper pair tunneling) lead 
to contributions of the form 

-^tt ~ ~ -^ctt sin(Ay -|- A(f) (96o) 

hi = - W sin( Ay - A(p) (966) 

with Ay = y 2 — yi the difference of the superconducting phases of the singlet pair 
potential on either side of the Josephson junction. The equilibrium configuration for 


zero current amounts tqH 

Ayeq = TT — arctan 


Icn - ^cu 


tan A(p 


(97) 


-6ctt + ^cU 

Thus, the system can act as a phase battery. For the half-metal case {Icn = 0) this has 
been noted using Eilenberger equations in Refs. and using Usadel equations 

in Ref. |224| . Going beyond the tunneling limit, we find a non-zero equilibrium phase 
difference between the superconducting leads, which shows a jump as function of Atp, 
indicating a first order transition (two competing minima of the free energy located 
at different values of phase difference Ay), which is illustrated in Fig. 19 (c). Note 
that the current-phase relation fulfills the symmetry I{—Ax,A(p) = —/(Ay, — Acp) 
following from the behavior of the current density under time reversal. The equilibrium 
configuration carries a non-zero spin current density Ig = 2I^(Axeq, Aip) ^ 0. The 
presence of a Josephson current at zero superconducting phase difference Ay was 
confirmed in calculations using Bogoliubov-de Gennes equations in a tri-layer geometry 
]. A study of a superconductor/ferromagnetic-insulator/superconductor junction on 


the surface of a three-dimensional topological insulator revealed similar effects |323] . 

The first-order processes described by Eqs. (96a)-(966) compete with a second order 
process, where two (or an even number of) pairs are transmitted simultaneously. In 
analogy to the so-called crossed Andreev reflection process, we call this process crossed 
pair transmission (see figure 20). For transmission of all pairs in equal direction, this 
process involves phases which are multiples of (Ay -|- A^p) -\- (Ay — Ap) — 2Ay, for 
which the dependence on Aip drops out. Consequently, the corresponding contribution 
to the Josephson current, Jcp, is independent of A^p and 7r-periodic. More generally, 
a transmission process of m spin-tt pairs and n spin-j-j- pairs involves a phase 
{m + n) Ay — n)A(p +{m + n)7r. Thus, the corresponding Josephson charge current 

can be written as 

2e + n){-iy 


I = 


\n+m 


Imn sin[(m -I- n) Ay + {m - n)A(p] 


n)Aip\ 

(98 o) 

n)Aip\. 

(986) 


whereas the spin current is 

Is = h^ ^(m - n)(-l)”+"^Jm„sin[(m -F n)Ay -\- {m - n)A(p\ 

run 

As I-m-n = Imn, the factor t can be omitted if we restrict the summation to m > 0, and 
for m = 0 to positive n. Note that the spin-current is even present in the case Ay = 0, 

^^For zero phase difference a spontaneous current /(Ay = 0) = — {/ctt “ Ttt) sin(Ay) appears. 
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Figure 20. Schematic drawing of the crossed pair transmission process. Cooper 
pairs with opposite spin projection on the magnetization vector do not contribute to 
the spin supercurrent in strongly spin-polarized ferromagnets (top panel). However, 
due to magnetic inhomogeneity across the interfaces singlet-triplet conversion processes 
generate out of two singlet pairs two equal-spin triplet pairs with opposite spin 
polarization. These two equal-spin pairs are transmitted simultaneously in a crossed 
pair transmission process (bottom panel). Their phases correspond to Ay ± Aip. A 
geometric phase shift 2Aip appears between the two types of Cooper pairs. Copyright 
(2009) by the American Physical Society |172| . 


provided that A(p ^ 0. An example of such a case was discussed in ma, where one of 
the superconductors was replaced by an insulator, and a pure spin supercurrent remains, 
giving rise to a spin-Josephson effect. In this case only terms with m+n = 0 remain, with 
a pure spin supercurrents Ig = On the other hand, crossed 

pair transmission processes in equal direction have zero spin current and correspond to 
m = n, leading to a contribution to charge current of the form 2e Xm sin(2mAx). 

For illustrative purposes we consider the example of only first order terms and leading 
order crossed pair transmission terms, i.e. only terms associated with /oi,/io, hi, and 
h-i- Then 

//2e = — ho sin(Ax + Aip) — /qi sin(Ax — A(p) + 2hi sin(2Ax) (99a) 
/s/h = — ho sin(Ax + Aip) + /qi sin(Ax — A(p) + 2h,-i sin(2A(y9). (996) 

The last term in each of these equations corresponds to the crossed pair transmission 
process (in equal and opposite direction, correspondingly). The critical current density 
for positive and negative current bias differs for this case, similar as shown in Fig. (b). 
Furthermore, for sufficiently large hi multiple minima of the free energy as function of 
Ax appear, leading to the characteristic jump at a certain value of Aip as illustrated 
in Fig. (c). If A(p were continuously varied, a hysteresis would occur, typical for a 
first order phase transition. Finally, when one of the superconductors is replaced by an 
insulating material then 1 = 0 and the spin-Josephson current is h = 2hh -i sin(2A(p) 

m- 
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6.4- Recent experimental observations 

6.4- 1- Triplet supercurrents in half-metallic ferromagnets Half-metallic ferromagnets 
are especially promising for applications, as they are fully spin polarized and thus give 
the largest spin hltering effect possible. The recent developments and applications for 
half-metallic ferromagnets have been reviewed e.g. in Refs. [3241132511326] . 

First indications for an coupling of superconductors through a half-metallic 
ferromagnet came from experiments by Pena et al. on ferromagnetic Lao.rCao.sMnOs, 
coupled with the high-temperature d-wave singlet superconductor YBa 2 Cu 307_5 in 
trilayers and superlattices |327[ I328j . These authors found a long-range proximity effect 
over a length scale of 100 nm. 

In 2006, Keizer et al. reported a triplet supercurrent in the half metal Cr 02 |329) . 
The authors used NbTiN as superconducting electrodes and Cr 02 grown on Ti 02 . 
After it had proved difficult to reproduce the effect, it was in 2010 when finally Anwar 
et al. in the group of Aarts reported that they had found long-range supercurrents 
through Cr 02 grown on AI 2 O 3 (sapphire) |330) . The current was observed over a 
distance of 700 nm between two superconducting amorphous Mo 7 oGe 3 o electrodes. The 
effect was interpreted in terms of odd-frequency pairing correlations (although really 
the experiment only indicates equal-spin pairing correlations and does not provide 
experimental insight about the orbital and frequency symmetry of the pairs). However, 
the group was not able to find the long-range supercurrent in Cr 02 grown on Ti 02 , 
the material used by Keizer et al. They noted that their Cr 02 films showed a uniaxial 
anisotropy in contrast to the biaxial anisotropy present in the Keizer et al. samples 

[ssniiHsi]. 

There had been criticism that based on previous point contact spectroscopy 
experiments, Cr 02 might not be fully spin-polarized. These data were summarized 
by Lbwfander et al. in 2010, and it was shown that all point contact data on Cr 02 
are indeed consistent with full spin polarization, if the theory is extended to take into 
account realistic interfaces lEa. 

Anwar et al. have in a subsequent study extended their work |332] and find that 
if a Ni/Cu sandwich between the Cr 02 film and the Mo 7 oGe 3 o electrodes is used, a 
long-range supercurrent is observed also for Ti 02 substrate over a distance of almost 
a micrometer. The critical current density is 100 times larger in these Ni/Cu/Gr 02 
junctions on Ti 02 substrate than if sapphire is used as a substrate, and is comparable 
to that observed by Keizer et al. This proves that in addition to the substrate on which 
Gr 02 is grown, the interface characteristics between Cr 02 and the superconductor play a 
decisive role, and spin-mixing due to misaligned spins as predicted before by theory |155] 
is the crucial ingredient for explaining the singlet-triplet conversion in these structures. 

Further evidence for triplet supercurrents came from a report by Sprungman et al ., 
who found a Josephson effect with ferromagnetic Cu 2 MnAl Heusler barriers |333) . It was 
observed that the critical current density versus temperature shows a pronounced peak 
around 4.5 K (electrodes were Nb), decreasing towards lower and higher temperatures. 
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Figure 21. (a) Schematic diagram of the Josephson junctions used in the work 

of Khaire et al.^ shown in cross-section. From Klose et al. |337| . Copyright (2012) 
by the American Physical Society, (b) Product of critical current times normal state 
resistance, /ci?N, as a function of total Co thickness, Dqq = 2dco- Red circles represent 
junctions with F’=PdNi and dpdNi =4 nm, whereas black squares represent junctions 
with no F’ layer (taken from Ref. |25Q) b As Dqo increases above 12 nm, IqR^ hardly 
drops in samples with PdNi, but drops very rapidly in samples without. (The solid 
line is a fit of the data without PdNi to a decaying exponential, from Ref. [250].) From 
Khaire et al. [M]. Copyright (2010) by the American Physical Society. 


Such a maximum was predicted in |155[ I133[ I215j . However, there is some intrinsic 
gradient of the degree of L2i-type Heusler structure ordering inside the Heusler layers, 
with a low degree if order at the interfaces and higher degree of order in the interior 
|333) ■ For that reason, it is too early to decide about the origin of this maximum in 
/c(T) in these experiments. 

Another half-metallic ferromagnet is LaoyCao.sMnaO, which has been studied 
in the context of YBa 2 Cu 307 /Laa;Cai_a;Mn 03 multilayers by Pena et al. in 2005. 
Kalcheim et al. |33^ performed scanning tunneling spectroscopy experiments on 
Lao. 7 Cao. 3 Mn 30 film epitaxially grown on superconducting Pri, 85 Ceo.i 5 Cu 04 , and find 
long-range penetration of superconductivity into half-metallic LaoyCao.sMnaO. Visani 
et al. |335] find quasiparticle and electron interference effects in the conductance across 
a Lao, 7 Cao, 3 Mn 30 /YBa 2 Cu 307 interface that demonstrate long-range propagation of 
superconducting correlations across the half metal. The effect is interpreted in terms 
of equal-spin Andreev reflections (equal-spin Andreev reflection was introduced in Ref. 
[T65] under the term “spin-flip Andreev reflection”, or SAR). The peculiarities of Andreev 
reflection at a half-metal interface with a singlet superconductor are connected with 
the spin-active interfaces that unavoidably will be involved, and have been clarified in 
Appendix C of Ref. |215j . and in Refs. |165[ 11581116711168] . 

6.4-2. Multilayer converter Experiments utilizing a multilayer geometry, in some sense 
manufacturing an inhomogeneous magnetization profile by using different materials for 
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the various layers, led to a successful generation of triplet supercurrents in a thick 
Co layer by the group of Birge in 2010 lam. In these experiments two layers of 
ferromagnetic alloys were used to generate triplet amplitudes, which then are passed as 
long-range supercurrent through a strongly spin-polarized ferromagnet. The crucial role 
of the misalignment of the magnetizations was demonstrated. With this breakthrough 
a simple and reliable way was found to produce long-range triplet supercurrents. The 
geometry is shown in Figure [^together with the results. It makes use of the previously 
acquired knowledge of the group that a Co/Ru/Co trilayer instead of a single Co layer 
will give much better control over the junction behavior [see Figure [^b)] |25()j . When 
the extra outer layers are present, a long-range effect persists, whereas if they are 
absent, there are only short-range supercurrents present |3()ll 1336] . It was later shown 
experimentally that the triplet supercurrent is enhanced up to 20 times after the samples 
are subject to a large in-plane magnetizing field aligning the two Co layers perpendicular 
to the magnetizations of the two thin outer layers [HHZ], exactly as predicted e.g. 
in [1551 11071 1296] . In Ref. |338j a triplet supercurrent in S/F’/F/F’/S Josephson 
junctions was studied, where S was superconducting Nb, F’ a thin Ni layer with in¬ 
plane magnetization, and F a Ni/[Co/Ni]„ multilayer with out-of-plane magnetization. 
It was found that the supercurrent decays very slowly with F-layer thickness and is much 
larger when the F’ layer are present than if not. This confirmed that the spin-triplet 
supercurrent is maximized by the orthogonality of the magnetizations in the F and F’ 
layers. 

Theoretical treatments of the geometry used by Khaire et al. have appeared in 
Refs. |300l I302[l283j . It was subsequently demonstrated experimentally that the critical 
current scales linearly with area in magnetized junctions, confirming the homogeneity 
of the superconducting phase difference across the junction area |339j . 

6.4.3. Long-range effects in ferromagnetic nanowires There have been also reports 
for a long-range proximity effect in long ferromagnetic nanowires. An early overview 
over studies of ferromagnetic nanowires with superconducting electrodes was given by 
Petrashov et al. |34()j . The controlled growth of nanowires has improved over the 
last decade, allowing to create crystalline nanowires of several hundreds nanometers 
length. Wang and coworkers used a single-crystalline ferromagnetic Co nanowire and 
reported zero resistance in wires up to 600 nm length m- A theory for a long- 
range singlet proximity effect in ferromagnetic nanowires was given by Konschelle et 
al. |342j . Takei et al. developed a theory of a ferromagnetic nanowire-superconductor 
proximity structure based on Rashba spin-orbit coupling in the barrier that induces 
p-wave superconductivity in the ferromagnet [MH] . 

Almog et al. |344j report measurements of the dynamical conductance of a double 
In/Co/In device, where two Co wires are connected to superconducting In electrodes, 
running parallel less than a coherence length apart from each other. They find a 
spin polarization of the superconducting order parameters at the interfaces which 
depends on the relative spin polarization of the two wires. Colei et al. |345j study 
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a similar structure with two superconducting electrodes being bridged by two parallel 
ferromagnetic wires forming an SFFS junction. They concentrate on the phenomenon 
of crossed Andreev reflection. At low temperatures and excitation energies below 
the superconducting gap, they find that the resistance corresponding to antiparallel 
alignment of the magnetization of the ferromagnetic wires is higher than that of parallel 
alignment. Spin-dependent interface scattering was found to be important in order to 
understand these findings [346] . 


6 .4-.4-. Spiral magnetic order A realization of the ideas by Bergeret, Volkov, and 
Efetov was reported in Ref. P7| . by employing the conical ferromagnet holmium. 
A long-range proximity effect was observed. In this work, superconducting phase- 
periodic conductance oscillations in ferromagnetic Ho wires in contact with conventional 
superconductors were measured. The distance between the interfaces was much larger 
than the singlet superconducting penetration depth. 

A triplet supercurrent has been found by Robinson, Witt, and Blamire with a 
setup using a cobalt layer sandwiched between two holmium layers, constituting another 
breakthrough in 2010 [282] . The main results of this work are shown in Fig. 22, Without 
the Ho layers between the Nb and Co (or with the Ho replaced by Rh) the characteristic 
voltage IcRn shows fast oscillations and a short-range decay over six orders of magnitude 
on a scale of 10 nm. If additional Ho layers are inserted between Nb and Co, the 
characteristic voltage decays very slowly, by about an order of magnitude on a scale of 
50 nm. The effect is spectacular especially for Co thicknesses above 10 nm. Theoretical 
treatments of this geometry have been given in Refs. |349[ 1350] 1185] . 

Further reports concentrate on creating artificial non-collinear magnetic structures 
that can be brought in contact with superconductors. In Ref. ESB it is reported that a 
new structure using an exchange-spring magnet was fabricated that can be tuned from 
a collinear to a non-collinear state by a rotating external magnetic field in a controllable 
way. With this setup the authors found an increase of superconductivity in the structure 
(transition temperature and conductance) when going to a non-collinear state. 

Robinson et al. studied Nb/Fe/Cr/Fe/Nb junctions where the thickness of the 
Cr layer determines the relative alignment of the Fe layers, and find a substantial 
enhancement of the critical Josephson current when a non-parallel configuration was 
realized [352] . 

Witt et al. find that in structures with a conical magnetic spacer layer of Ho the 
critical Josephson current decays exponentially with layer thickness, corresponding to a 
coherence length of 4.34 nm in Ho, and without showing any oscillatory behavior [353] . 
Comparing this with the mean free path of 0.28-0.87 nm, they conclude that their Ho 
structures are in the dirty limit. 

A realization of field-tunable in-plane Bloch domain walls in the rare-earth magnet 
gadolinium if grown between non-collinearly aligned ferromagnets was suggested in 
Ref. [354] . It was found that supercurrents flow through magnetic Ni/Gd/Ni nanopillars, 
the magnitude of which strongly depends on the domain wall state in Gd. The authors 
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Figure 22. (A) Slow decay at 4.2 K in the characteristic voltage IcRn of 

Nb/Ho(4.5 nm)/Co( (ico)/Ho(4.5 nm)/Nb junctions (blue circles) and a Nb/Ho(10 
nni)/Co(dco)/Ho(10 nni)/Nb junction (green circle) versus Co barrier thickness 
(dco)- Inset: Comparative data (black circles) from |348| showing the behavior of 
Nb/Rh/Co/Rh/Nb junctions. The oscillating curves in the inset and main panel are 
theoretical fits to the experimental data in the inset, as described in [348]. (B) The 
conical magnetic configuration of idealized Ho below its Curie temperature (20 K), 
showing an antiferromagnetic spiral rotating in-plane 0 = 30° per atomic plane and 
pitched a = 80° out-of-plane. The moments (arrows) rotate about the surface of 
a cone with the spiral wavelength, A, corresponding to a Ho thickness of 3.4 nm. 
(C) Device layout consisting of two superconducting Nb electrodes coupled via a 
Ho/Co/Ho trilayer. From Robinson et al. |282) . Reprinted with permission from 
AAAS. 


explain this result in terms of the inter-conversion of triplet and singlet pairs, the 
efficiency of which depends on the magnetic helicity of the structure. 

For a recent review by Blamire and Robinson see [355] . 

7. Modern developments 

In the following I selectively give examples for exciting modern developments of the 
field, without claiming to cover the entire spectrum of activities. 

7.1. Spin-valve devices and spin-filter junctions 

Although currently the research still concentrates on studying fundamental questions, 
applications can be imagined in various ways. The most obvious application would be 
the development of a spin-valve device. Re-entrant phenomena are predicted not only 
as function of thickness of the adjacent ferromagnetic layers, but also as function of 
other parameters, e.g. the degree of magnetic inhomogeneity, as in ferromagnets with 
spiral order. A high degree of experimental control has been achieved in producing 
large batches of samples with gradually varying thicknesses, allowing for a tailoring of 
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the Fulde-Ferrell-Larkin-Ovchinnikov state |656) . A controllable Josephson spin-valve 
has been realised recently ES7|. 

Superconducting spin-filter tunnel junctions have been studied recently experimen¬ 
tally |358) and theoretically |359) . In Ref. |358) , GdN barriers are used as ferromagnetic 
insulator barriers, and it is shown that the field and temperature dependence of the crit¬ 
ical Josephson current is strongly modified by the ferromagnetic insulator. It is found 
that the strong suppression of Cooper pair tunneling by the spin filtering of the barrier 
can be modified by magnetic inhomogeneity in the barrier. In theoretical work [359] it 
was also demonstrated that the differential conductance may exhibit peaks at different 
values of the voltage depending on the polarization of the spin filter, and the relative 
angle between the exchange fields and the magnetization of the barrier. 

1.2. Flux Qubits and semifluxons 

Superconducting electronics is undergoing a revival at present, with pressing need for 
cryogenic memory for flux quantum logic applications as well as new superconducting 
Qubit applications. Metal spintronics is already widely employed in computer hard 
disc technology, and superconducting spintronics may lead to energy efficient memory 
and logic for supercomputing applications. Superconducting digital single-flux-quantum 
circuits using superconductor-ferromagnet-superconductor sandwich technology to 
insert tt-J osephson junctions into a circuit is promising due to its high operation speed 
and low energy consumption [360] . Such systems also solve the problem of high element 
densities on-chip for operating and holding magnetic flux quanta. A combination of 
magnetic Josephson junctions and conventional Josephson junctions can be used to form 
addressable memory cells, energy-efficient memory periphery circuits and programmable 
logic elements [36111362] . 

Another interesting development is semifluxon physics in zero-7r Josephson 
junctions. Fluxons are traveling, solitonic waves of magnetic flux in long Josephson 
junctions, created by external magnetic fields. In a zero-7r-junction a Josephson vortex 
of fractional magnetic flux is pinned at the zero-yr-boundary. A memory cell based on 
a (p-Josephson junction has been suggested [363] . where writing is done by applying a 
magnetic field and reading by applying a bias current. Storage at low temperatures is 
passive, without any bias or magnetic field applied. 

A high-frequency cryogenic generator operating at about 200 GHz, and based 
on flipping a semifluxon in a Josephson junction has recently been demonstrated in 
Ref. [366] . A tt junction is artificially created by current injection using a technique by 
Ustinov [364] . giving rise to a semifluxon [365] . A conversion efficiency of ~ 10% of dc 
input power to ac output power was achieved, including the losses on the way from the 
generator to the on-chip detector. This type of Josephson oscillator is comparable with 
those based on flux-flow, with advantages like small size and insensitivity to injection 
current. 

Realization of a O-tt junction in an atomic bosonic quantum gas has been proposed 
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in Refs. |667[ 1668) . In Ref. |668) it is suggested that two-state atoms in a double-well 
trap are coupled and an all-optical O-vr Josephson junction is created by the phase of a 
complex-valued Rabi frequency, exhibiting modes similar to semifluxons. It is suggested 
that pairs of semifluxons can be created by starting from a flat-phase state in long, 
optical O-TT-O Josephson junctions formed with internal electronic states of atomic Bose- 
Einstein condensates |369) . 

7.3. Non-equilibrium quasiparticle distribution 

A particular exciting subject is the possibility to utilize non-equilibrium quasiparticle 
distribution in conjunction with quantum coherence in so-called Andreev interferometer 
geometries. These devices were introduced by Petrashov, Antonov, Delsing, and Claeson 
in 1994 |67()) . A superconducting wire (e.g. aluminum) is attached at two points to a 
normal mesoscopic conductor (e.g silver or antimony) to form a loop. The conductance 
of the wire then oscillates as function of the magnetic flux through the loop. The 
oscillations are attributed to phase transfer from the superconducting condensate to 
normal electrons via Andreev reflections at the N-S interfaces. The so-called tt-SQUID 
employs an idea of Volkov IHTH, where by applying a control voltage directly to a 
normal metal the electron distribution function in the normal metal is modified. This 
non-equilibrium distribution spreads to the superconductor and influences the transport 
in a nonlocal way. This allows for a direct control of the supercurrent through the 
device, including switching between zero-junction and 7r-junction behavior. Such a 
device was realized experimentally by Morpurgo, Baselmans, van Wees, and Klapwijk 
in 1998-99 |372i 1373) . Corresponding effects are intricately connected with non-locality. 
Cadden-Zimansky, Wei, and Chandrasekhar combined an Andreev interferometer 
arrangement with injection of non-equilibrium excitations into a normal wire, giving 
rise to a current that adds to the supercurrent in the Andreev interferometer and leads 
to a voltage at the contacts. Due to current conservation each current influences the 
other, and together they allow for tuning of the contact voltage by changing the flux 
through the Andreev interferometer that controls the supercurrent Bia. Coherent 
voltage oscillations as functions of external flux are the result. 

The combination of non-equilibrium quasiparticle distribution with the Josephson 
effect seems particularly exciting, and is largely unexplored so far. Recent studies by 
Bobkova and Bobkov |376[ 1377) address the problem of spin-dependent non-equilibrium 
quasiparticle distribution. It is found that the interplay between the spin-dependent 
quasiparticle distribution and the triplet superconducting correlations induced by the 
proximity effect between the superconducting leads and ferromagnetic elements of the 
interlayer leads to the appearance of an additional contribution to the Josephson current. 
The interplay between short-range and long-range proximity effect is elucidated, and a 
long-range penetration of opposite-spin Cooper pairs under non-equilibrium conditions 
is proposed. An increase of the critical Josephson current by few orders of magnitude 
as a result of the non-equilibrium population in the ferromagnetic layer is suggested. 
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Long-range spin and charge accumulation in mesoscopic superconductors with Zeeman 
splitting was recently also studied by Silaev et al. [HZH]. 

A quantum interference transistor involving textured ferromagnets was suggested 
in Ref. |379) . It was shown that such a device acts as an ultra-sensitive magnetometer 
and allows for singlet-triplet switching by tuning a bias voltage. A combination of two 
spin injectors with an SFS Josephson junction was discussed in Refs. |,H8()[ 1881) . 


7.4- Dynamical effects 


Dynamical effects in superconductor-ferromagnet structures have moved in the focus of 
interest recently. In particular, it is clear that magnetization dynamics will be linked 
closely with Josephson dynamics, leading to potentially new effects. The geometric 
phases discussed in section |6.3.4[ when time dependent, are expected to lead to spin 
accumulation effects via an imbalance between the equal-spin pair electrochemical 
potentials, very similar to the voltage linked to a dynamical superconducting phase 
in the ac Josephson effect nza. 

The Josephson current in a diffusive superconductor-ferromagnet-superconductor 
junction with processing bulk magnetization was calculated in Ref. |882) . It was found 
that when the junction is phase biased, a dc Josephson current without ac component 
can still flow under this non-equilibrium condition. Long-range triplet amplitudes are 
induced by the processing magnetization. 

In Refs. |383[ 1384] non-equilibrium effects in a Josephson junction with two s- 
wave singlet superconducting leads coupled via a processing spin in a quantum dot 
was examined. An external magnetic field leads to a Larmor precession of the spin, 
rendering the magnetically active interface time-dependent. The authors analyze the 
non-equilibrium population of Andreev sidebands and dynamical spin currents. It is 
found that the supercurrent is enhanced and the critical Josephson current density shows 
a non-monotonous behavior as function of temperature, accompanied by a corresponding 
change in spin-transfer torques acting on the precessing spin. 

Supercurrent-induced magnetization dynamics in superconductor-ferromagnet 
Josephson junctions was explored in Refs. [38511386) 1387) 138811389) 1390) . It is found 
that the spin supercurrent can induce magnetization switching that is controlled by the 
superconducting phase difference. The authors in Refs. |388[ 1389) confirm the finding 
of Ref. dza, that the effect of chiral spin symmetry breaking of the structure leads to 
additional geometric phases that allow for the stabilization of a cp-junction. 

In Ref. m the dynamics of superconductor-ferromagnet-insulator-ferromagnet- 
superconductor (SFIFS) junctions with a thin ferromagnetic layers investigated. The 
coupled dynamics of the magnetization and the Josephson phase leads to Josephson 
plasma waves coupled to oscillations of the magnetization, affecting the form of the 
current-voltage characteristics in weak magnetic fields. 
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7.5. Spin-orbit coupling and topological materials 


Spin-orbit effects play an important role whenever the inversion symmetry is broken, 
either in the bulk material (when it is lacking a center of inversion), or at interfaces and 
surfaces |392j . A prominent example for spin-orbit effects at surfaces is the Rashba- 
Bychkov spin-orbit coupling |393] . and in the bulk an example is the Dresselhaus 
coupling due to bulk-inversion asymmetry |394j . 

Whereas in materials with a center of inversion all band-diagonal matrix elements 
of the spin-orbit coupling vanish, this is not the case in non-centrosymmetric materials. 
Spin-orbit coupling in non-centrosymmetric materials is strongly enhanced due to band- 
diagonal contributions, leading to a splitting of the Fermi surface into spin-orbit bands, 
sometimes also called ‘helicity bands’ (although a well defined helicity is only in special 
cases present). The kinetic part of the Hamiltonian in a one-band model has the form 


"Hkin = 


E E 


E(k) + g(k) ■ 


( 100 ) 


with the spin-orbit vector being odd in k: g(—k) = — g(k). In figure 23 various 
spin-orbit vector fields compatible with the point group of the crystal are visualized 
on a hypothetical spherical Fermi surface. All three cases are relevant for materials: 

for CePtaSi |395| . CeRhSia |396j and CeIrSia |397j . for Y 2 Ca [398] . and O for 
Li 2 (Pdi_a;Pta;) 3 B [399] . The kinetic part of the Hamiltonian can be diagonalized, leading 
to the above-mentioned helicity bands. The spin-orbit interaction locks the orientation 
of the quasiparticle spin with respect to its momentum in each band. 

One interesting aspect of a non-centrosymmetric ferromagnetic Josephson junction 
is the modification of the Josephson relation from an odd function in the snperconducting 
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Figure 23. Various types of spin-orbit coupling in crystals, relevant for topological 
superconductivity. The crystal point group sets restrictions to the allowed spin-orbit 
vector fields. The lowest order expansion terms in crystal momentum are shown. 
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phase difference to a current-phase relation where this symmetry is broken. For example, 
near the critical temperature the current-phase relation I = Ic sin(Ay) can be modified 
to 

/ =/cSin(Ax - (po), (101) 

where Ay is the superconducting phase difference, and (po is a phase shift proportional 
to the magnetic moment perpendicular to the gradient of the asymmetric spin-orbit 
potential Hm. The possibility of a (po-j unction has been already taken into account by 
Josephson ISH, and was later considered in Josephson junctions involving unconventional 
superconductors [1001. Examples in Josephson junctions with half-metals and strongly 
spin-polarized ferromagnets we discussed in Section |6.3.4[ A similar effect has also been 
predicted for Josephson junctions with a spin-polarized quantum point contact in a two- 
dimensional electron gas with spin-orbit coupling in an external magnetic field |402j . 
The spin-dynamics with such a (po Josephson junction has been discussed in Ref. |3Hb) . 
The competition between a Zeeman interaction and a Rashba spin-orbit interaction 
has been in the center of attention since a while. In connection with one-dimensional 
quantum wires, an anomalous cpo shift was predicted in Refs. |4()31I4()4) . In a model for a 
mesoscopic multilevel quantum dot the conditions for an anomalous Josephson current 
Eq. (101) were found to be a finite spin-orbit coupling, a suitably oriented Zeeman 
field, and the dot being a chiral conductor [40511406] . In Refs. [32211407] an anomalous 
Josephson current was predicted in junctions coupled with a two-dimensional electron 
gas exhibiting coexistence of spin-orbit coupling and Zeeman field. 

A recent development concerns a gauge-covariant approach to establish the 
transport equations, treating the charge and spin degrees of freedom on equal footing. 
In this approach both the electromagnetic and spin interactions are described in terms 
of U(l) Maxwell and SU(2) Yang-Mills equations, respectively |4()H[ 1409] . The starting 
point is a quasiclassical expansion of the microscopic Gor’kov equations in a gauge- 
covariant manner jnniiiii]. The idea is that one can re-write a Hamiltonian of the 
type 


as 


TJkin — 


TJkin — 


Px+Py 

2m 


- /i + h • (T Us.o.(z X P) • cr 


iPx iPy T 


+ 




mv: 


( 102 ) 


Th-cr, (103) 


2m 2m 2 

which motivates to introduce = —mVs.o.<yy/h cmd Ay = mVs.o.crx/^ as two non-abelian 

gauge potentials. In this approach, the covariant derivative 

= dRG + i[A,G] (104) 

with the gauge potential A is associated with the gauge field 


dyAy djjAy 


f [Ay,t, A^, 


(105) 


which is antisymmetric in its indices, and has only one component, Yxy = 
—2(mVs.o.Yo'z/■ This leads in the transport equation to the modification 

ihvf ■ Vng ifwf • Vr^ -F imvl^ d^ {<yz-,9} 


(106) 
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with the derivative along the Fermi surface (which is assumed a circle for simplicity). 
Using this modified equation, the role of spin-orbit coupling as source of long-range 
triplet proximity effect in superconductor-ferromagnet structures has been explored in 
Ref. m- Furthermore, the connection to ipo Josephson junction behavior has been 
elucidated [412[I413| . It should be, however, cautened that the modification in Eq. (106) 
does not contain all terms of the same order as the additional term. A quasiclassical 
theory of disordered Rashba superconductors has also been proposed in Ref. Em- 

One way to manipulate spin in spintronics devices is the so-called spin Hanle effect 
prai?rai4T7| . which describes the coherent rotation of a spin in an external magnetic 
field. A recent theoretical study of this effect is given in Ref. |418j . It is demonstrated 
that superconductivity can strongly influence the coherent spin rotation, depending on 
the type of spin relaxation mechanism being dominated either by spin-orbit coupling or 
spin-flip scattering at impurities. 

Another manifestation of coherence in systems with spin-orbit coupling is the 
Aharonov-Casher effect Em, leading to a phase on the Josephson current through a 
semiconducting ring attached to superconducting leads |42()) . This effect is the charge- 


spin dual effect to the Aharonov-Bohm effect Em, which describes the relative phase 
shift between two charged particle paths enclosing a magnetic flux. Both effects are 
manifestations of geometric phases acquired by the quantum mechanic wave function 
under adiabatic changes 112211122]- The particle’s spin acquires such a geometric phase in 
systems with spin-orbit interaction |424[ 1425114261142711428L142911430) . The Aharonov- 
Casher effect was, e.g., observed experimentally in ring structures of HgTe/HgCdTe 
quantum wells Em- In superconducting Josephson rings, or closed Josephson junction 
arrays, such phases lead to oscillations of the Josephson current due to the Aharonov- 
Casher phase [4321 14201 14331 1434) . The effect allows for control of the Josephson 
current through the control of the Aharonov-Casher phase by the gate voltage. Thus, 
this effect is a promising candidate for realizing new types of controllable devices in 
superconducting spintronics based on geometric phases. 

Recently Mironov et al. BSl studied double path interference during Cooper pair 
transport through a single nanowire with two conductive channels. It is found that multi¬ 
period magnetic oscillations appear due to quantum mechanical interference between 
channels affected by spin-orbit coupling and Zeeman coupling. The model is relevant to 
recent observations of interference phenomena in Bi nanowires 


A Josephson junction containing a spacer with strong spin-orbit interaction was 
considered in Ref. P7) . A nonlinear dynamical coupling between magnetic moment and 
charge current was found, and magnetic torque and charge pumping was investigated 
in such a system. The intricate coupling between spin and charge currents in 
systems with strong spin-orbit coupling is previously known from non-centrosymmetric 
superconductors, where spin-polarized Andreev states play a prominent role 1138]. 

This connects to the current hot topic of topological materials, in particular 
topological insulators and superconductors, both systems with strong spin-orbit 
coupling. As an example for the extremely rich plethora of effects involving topological 
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materials we mention here the possibility to observe chiral Majorana modes in one¬ 
dimensional channels built at a superconductor/ferromagnetic-insulator/superconductor 
junction on top of a topological insulator |323j . where a ipQ junction is predicted. The 
effect is interpreted as tunneling process between two Majorana edge channels at the 
two interfaces between the superconductor and the ferromagnetic insulator. Majorana 
fermions have the property of being their own antiparticle, i.e. their field operators fulfill 
7a = 7a |439[ 1440] . A non-trivial superconducting phase is also obtained in proximity 
junctions involving semiconductors with Rashba spin-orbit coupling and a time-reversal 
symmetry breaking Zeeman term in the Hamiltonian |441l 1442] . 

The combined effect of spin-orbit interaction, magnetic field, and Coulomb charging 
for a multilevel quantum dot tunnel contacted by two superconductors was analyzed in 
Ref. [443] . Majorana bound states in a double dot variant of this system are predicted to 
leave a clear signature in the 27r-periodic current-phase relation. Majoranas in spin-orbit 
coupled ferromagnetic Josephson junctions were investigated in Ref. |m] , were it was 
shown that two delocalized Majorana fermions with no excitation gap appear in a tt- 
junction. Josephson currents through Majorana bound states in topological insulators 
have been studied in Refs. |445[ 14461 1447] 14481 1449] 1450] 14511 1452] 14531 1454] . Such 
Josephson junctions carry Jvr-periodic bound states |445j . It has been shown that under 
certain conditions this periodicity manifests itself by an even-odd effect in Shapiro steps 
[4471 1453] . In addition, a peak in the current noise spectrum at half the Josephson 
frequency has been predicted |454j . 

The diverse spectrum of effects and phenomena in hybrid systems between singlet 
superconductors and topological insulators (see e.g. Refs. [455114561145711458] 1459111371 
146011461] ) are promising examples of how the field can be brought forward, playing a 
prominent role in various modern developments at the forefront of international research. 
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